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AESTRACT 


Hvdrodynamical  theories  on  meanders  of  the  Gulf  Stream  and 
similar  western  boundary  currents  are  critically  reviewed.  In  Section  1, 
observed  facts  of  meanders  east  of  Cape  Hatteras  are  described  including 
data  obtained  in  Operation  Cabot  and  Gulf  Stream' 60.  In  Section  2,  long 
v/aves  along  a  shearing  horizontal  current  in  a  rotating  barotropic  system 
are  discussed,  based  on  Haurwitz  and  Panofsky's  model.  In  Section  3, 
stability  of  long  waves  in  a  uniform  current  of  a  two-layer  ocean  is 
treated  from  Stommel's  and  Iida  and  Ichiye's  theory.  In  Section  4,  theories 
on  waves  on  a  shearing  current  in  the  upper  layer  of  a  two-layer  including 
Duxbury,  Iwata  and  Yasui  are  presentedc  In  Section  5,  general  criteria 
on  stability  of  a  baroclinic  flow  are  discussed,  including  Pedlosky's 
application  of  Lin's  theorem.  In  Section  6,  the  Rossby  waves  in  a  baro¬ 
tropic  and  baroclinic  ocean  are  explained  from  various  authors  including 
Fofonoff,  ivloore,  II' in  and  Kamenkovich,  R  elander  and  Ichiye  and  insta¬ 
bility  of  waves  in  an  eastward  flow  is  discussed  as  a  cause  of  meanders.  In 
Section  7,  applications  of  Rossby' s  constant  absolute  trajectory  theory  n  to 
meanders  by  Saint-Guily,  Warren  and  Ichiye  are  presented  and  Greenspan's 
criterion  on  stability  of  a  western  boundary  current  is  derived.  In 
Section  8,  the  cross  current  momentum  eddy  transfer  determined  by 
Webster  and  Ichiye  for  the  Gulf  Stream  and  the  Kuroshio,  respectively 
and  the  eddy  heat  transfer  determined  by  Oort  for  the  Gulf  Stream  are 
discussed.  In  Section  9,  possibility  of  testing  some  theories  on  meanders 
by  use  of  a  rotating  tank  is  discussed  on  the  basis  of  Faller's  explanation 
of  similarity  between  p-effect  and  variable  depth.  In  the  appendix,  the 

observed  facts  of  meanders  of  the  Kuroshio  are  presented  from  results  of 
Masuzawa,  vloriyasu,  Uda  and  Ichiye. 
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1.  INTRODUCTION 


One  of  the  most  interesting  features  of  the  Gulf  Stream  is  its 
meandering  which  seems  to  resemble  in  some  respects  the  meanders  of  a  river 
in  its  flow  pattern.  Iselin  (1936)  noticed  that  the  horizontal  distribution  of  iso~ 
therms  and  dynamic  topography  of  the  western  North  Atlantic  indicate  wavy 
patterns  along  the  path  of  the  Gulf  Stream  to  the  south  of  Grand  Banks  after  the 
G'df  Stream  leaves  the  continental  shelf.  Church  (1937)  confirmed  that  the 
position  of  the  Gulf  Stream  between  Mont  auk  Point  and  Bermuda  is  highly  vari¬ 
able  from  surface  temperature  data  obtained  with  thermograph  installed  on 
liners  „ 

Observational  techniques  for  study  of  meanders  of  the  Gulf  Stream 
were  substantially  improved  by  introduction  of  BT,  GEK,  Sea  Sampler  and 
STD  Recorder  together  with  Loran  which  provides  the  necessary  accuracy  in 
navigation.  These  instruments  have  made  it  possible  to  measure  respectively 
the  temperature  in  the  upper  mixed  layer  and  surface  currents  on  board  a  ship 
moving  with  a  cruising  speed,  and  therefore,  to  collect  the  data  quasi- 
synopticallyc  Ford  and  Miller  (1952)  showed  that  the  Gulf  Stream  meander  - 
ings  developed  into  a  large  anti-cyclonic  loop  of  about  150  nautical  miles 
diameter  in  the  area  between  62°  to  66°  W  and  38°  to  41®  N. 

These  observations  were  taken  by  one  ship  on  the  edge  of  the  Stream 
for  1050  miles  between  Cape  Hatteras  and  63®  W  from  November  15th  to 
December  6th,  1948,  using  these  new  techniques  and  instruments.  Besides, 
they  could  show  that  within  a  period  of  two  weeks  the  mean  axis  of  a  150-mile 
segment  of  the  Stream  between  66®  and  68®  V’  changed  its  direction  by  90°, 


. 
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suggesting  that  the  propagation  speed  of  the  meander  crest  (northward  pro¬ 
trusion)  was  about  11  n.  miles  per  day. 

However,  the  first  extensive  survey  of  the  Gulf  Stream  deploying 
almost  simultaneously  six  research  vessels  equipped  with  these  instruments, 
was  made  between  Cape  Hatteras  and  the  Grand  Banks  during  June  6-23,  1950 
and  is  called  "Operation  Cabot"  (Fuglister  and  V/orthington,  1951).  During 
"Operation  Cabot",  almost  daily  synoptic  plots  of  a  mean  temperature  of  the 
upper  200-meter  layer  were  made  with  BT  data.  These  data  indicate  that  the 
crest  of  the  first  wave  near  72®  W  in  meanders  of  the  Gulf  Stream  moved 
toward  the  east  at  a  rate  of  approximately  11  nautical  miles  a  day,  as  shown 
in  Fig.  1,  Also,  in  this  operation  it  was  observed  that  a  large  cyclonic  eddy 
consisting  of  cold  water  was  formed  at  about  59~6l®  W  and  37-40®  N  when  a 
large  loop  in  the  Gulf  Stream  became  cut  off  and  separated  from  the  Stream 
in  about  nine  days.  The  cut-off  eddy  was  elongated  in  the  direction  of  west 
northwest  to  east  southeast  with  a  long  axis  of  about  180  n.  miles  and  a  short 
axis  of  about  70  n.  miles.  It  was  not  known  from  this  operation  how  long  it 
took  for  the  trough  to  develop  Into  a  narrow  necked  loop.  Also,  the  movement 
and  duration  of  the  eddy  after  its  cut-off  from  the  Stream  and  statistical  natures 
of  such  eddies  such  as  typical  patterns,  and  place  of  occurrence  were  left  un¬ 
answered  by  this  study. 

Fuglister  (1951,  1955)  proposed  a  different  interpretation  of 
horizontal  temperature  fields  of  the  Gulf  Stream  and  advanced  a  hypothesis  of 
multiple  currents  which  consist  of  a  series  of  overlapping  currents  and  rela¬ 
tively  weak  countercurrents  instead  of  one  continuous,  branching  current  as 
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heretofore  generally  accepted.  According  to  this  interpretation,  an  idea  of 
a  single  current  of  the  Gulf  Stream  with  the  mear.derings  which  increased 
their  amplitudes  to  the  east  of  65°  VI  is  abolished.  Instead,  it  is  interpreted 
that  several  currents  with  relatively  straight  courses  flow  side  by  side  sepa¬ 
rated  by  v/eak  countercurrents  and  correspondingly,  the  isotherms  run 
approximately  in  zonal  directions  with  successive  cold  and  warm  belts.  In 
between  these  two  extreme  interpretations,  there  is  a  type  of  analysis  which 
retains  an  idea  of  a  single  current  but  manages  to  avoid  violent  and  improb¬ 
able  meandering  by  introducing  eddies.  Fuglister  applied  these  types  of 
analysis  to  the  temperature  field  at  a  200  m  depth  of  the  Gulf  Stream  between 
75°  and  52°  W  observed  by  three  ships  in  August,  1953  and  of  the  Kuroshio 
between  140°  and  165®  E  observed  by  a  multiple  ship  survey  in  August,  1933, 
The  idea  of  multiple  currents  proposed  by  Fuglister  was  based 
on  temperature  data  alone  and  the  overlapping  segments  were  described  as 
"shingles"  from  the  appearance  of  the  corresponding  isotherms.  Von  Arx 
and  others  (1955)  found  that  the  shoreward  edge  of  the  Gulf  Stream  between 
Florida  and  70®  W  was  not  a  continuous  front  but  consisted  of  shingles  with 
lengths  from  100  to  500  km.  They  tracked  the  shingles  both  v/ith  a  ship  and 
an  aircraft  which  was  flying  at  an  altitude  of  1500  feet,  scanning  the  sea 
surface  with  a  radiation  thermometer.  These  shingles  were  easily  recog¬ 
nized  by  the  strong  temperature  gradient  across  them,  the  change  of  sea 
state  and  sea  color  and  the  accumulated  Sargassum  and  flotsam.  In  each 
segment  these  characteristics  became  fainter  as  observations  were  pursued 
downstream  and  eventually  disappeared,  but  always  another  frontal  segment 
was  encountered  at  about  90®  left.  Particularly,  visual  observations  from 
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the  aircraft  indicated  that  a  fine  shingle  structure  interrupted  the  frontal 

segment  at  an  interval  of  about  10  km  and  usually  was  left-handed.  They 

proposed  a  possible  interpretation  of  the  shingle  structure  as  an  effect  of 

the  diurnal  tide  of  the  Gulf  of  Mexico  which  supplies  the  additional  water 

transport  of  10  x  10^  m^/  sec  to  the  Straits  of  Florida  during  a  twelve-hour 

ebb  period  and  withholds  the  same  amount  during  a  flood  period,  since  the 

tidal  waves  enter  the  Gulf  through  the  Yucatan  Channel  and  leave  it  through 

the  Straits  of  Florida,  the  averaged  outflow  through  which  amounts  to 
6  3 

30  x  10°  m  /sec.  They  speculated  that  the  water  entering  the  Gulf  during 
a  flood  period  is  warmed  there  and  flows  out  into  the  Florida  Current,  form¬ 
ing  warm  blobs  in  the  warm  core  of  the  current  further  north  and  shingle 
segments  on  the  shoreward  edge.  In  fact,  the  core  of  the  Gulf  Stream  with 
characteristic  velocities  of  3  to  6  knots  transports  a  parcel  of  water  for  a 
distance  of  130  to  260  km.  during  a  24-hour  period,  corresponding  to  the 
lengths  of  shingle  segments  and  blobs.  However,  the  warm  blobs  which  are 
50  to  100  krn  long  are  commonly  located  with  distances  of  100  to  200  km. 
between  centers  in  the  Kuroshio  near  Japan,  although  there  is  no  counter¬ 
part  of  the  Gulf  of  Mexico.  (Ichiye,  1955).  It  is  yet  to  be  determined  whether 
there  is  the  shingle  structure  on  the  shoreward  edge  of  the  Kuroshio. 

It  is  considered  that  the  Gulf  Stream  does  not  show  large  meander- 
ings  to  the  south  of  Cape  Hatteras,  However,  26  serial  sections  made  at 
right  angles  to  the  mean  axis  of  the  Stream  off  Onslow  Bay  from  the  16th  of 
vlay  through  June  7,  1953  indicated  that  the  velocity  components  lateral  to 
the  mean  axis  of  the  Stream  could  change  within  a  matter  of  a  day  or  two  from 
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a  diffluent  to  a  confluent  pattern  (Von  Arx  and  others,  1955),  Serial  sections 
off  Onslow  Bay  were  again  made  by  a  Woods  Hole  group  from  May  20th 
through  June  18,  1958.  These  data  together  with  those  obtained  in  the  Straits 
of  Florida  between  Miami  and  Gun  Cay  by  a  University  of  Miami  group  from 
1952  through  1958  were  analyzed  to  determine  kinetic  energy  balance  across 
the  Stream  by  its  meanders.  (Webster,  1961  a,  b).  Serial  sections  of  a 
similar  kind  were  made  seventy  times  along  the  30°  N  parallel  off  Jacksonville 
during  October  and  November,  1961  and  twenty  times  off  Cape  Hatteras  dur¬ 
ing  October,  1962  by  another  Woods  Hole  group  and  these  results  were  ana¬ 
lyzed  by  Oort  (1964)  to  determine  heat  and  mass  transports  across  the  Stream. 
These  two  analyses  will  be  discussed  in  a  chapter  below. 

In  the  spring  of  I960  for  over  a  period  of  two  and  one-half  months 
from  the  2nd  of  April  to  the  15th  of  June,  an  area  from  Georges  Bank  off 
Cape  Cod  to  the  east  of  the  Grand  Banks  and  from  the  continental  shelf  south 
to  the  latitude  of  Bermuda  was  surveyed  by  three  oods  Hole  research  vessels 
(Fuglister,  1963).  These  surveys  included  deep  hydrographic  sections  along 
the  whole  length  of  meridians  in  the  area  with  spacings  of  two  degrees  of 
longitude  and  with  distances  between  stations  ranging  from  twenty  miles  on 
the  continental  slope  to  sixty  miles  south  of  the  Gulf  Stream,  tracking 
neutrally-buoyant  floats  at  depths  between  2000  and  4000  m,  and  transponding 
surface  buoys  and  measuring  the  thermal  structure  to  a  depth  of  150  m  with 
towed  thermistors.  It  was  unexpectedly  found  that  during  the  entire  two  and 
one-half  months  the  large  meanders  changed  position  by  less  than  the  width 
of  the  current  and  that  every  observed  change  in  the  position  of  the  current 
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could  be  accounted  fer  by  lateral  shifts  cf  me  Stream  with  speeds  less  than 

2.5  miles  per  day.  The  operation  of  tracing  the  position  of  maximum  current 
along  the  Stream  by  use  of  GEK  so  that  it  registers  no  component  normal 

to  such  a  course  indicated,  on  many  occasions,  that  the  path  of  maximum 
current  was  not  a  single  smooth  curve  but  had  shingle  structure.  The 
feature  of  the  Gulf  Stream  in  this  survey  showed  that  the  major  current  flowed 
in  a  direction  a  little  north  of  east  without  meandering  from  69°  to  62,5°  W, 
but  it  turned  abruptly  northward  and  formed  a  large  loop,  centered  around 

61.5  °W,  Subsequently,  the  current  headed  due  south  for  a  distance  of  over 
200  miles  along  60,  5®  vY  forming  a  "sock"  or  a  cyclonic  loop  along  the  east 
side  of  which  the  major  current  turned  to  the  north  again  from  37°  N,  60.  5®W 
to  41 0  5°  N,  53°  V7 ,  An  eddy  with  cold  water  mass  seemed  to  form  at  the  'toe' 
of  the  "sock"  or  the  southwesterly  end  at  the  beginning  of  the  survey,  moving 
slowly  to  the  north  along  an  anti-cyclonic  curve  and  one  month  later  this  eddy 
was  observed  with  its  center  at  36°  50'  N;  64°  30'  W,  elongating  from  west 
northwest  to  east  southeast  with  a  major  axis  of  200  miles  and  a  minor  axis 

of  fifty  miles.  To  the  east  of  the  "sock",  the  current  formed  minor  meander- 
ings  whose  ere  sts  and  troughs  were  located  at  41 , 5°  N, ,  57.5°  W;  38,5°  N, 
55°  W;  40°  N,  54°  W  and  39°  N,  52,5°  V>;,  respectively,  Fuglister  (1963) 
indicates  that  the  Gulf  Stream  meanders  are  not  a  series  of  waves  gradually 
increasing  in  amplitude  from  west  to  east,  but  rather  a  quasi-stationary 


pattern  with  an  abrupt  increase  of  the  amplitude  near  62°  W  from  the  data 


of  various  cruises  since  1947  in  an  area  between  60®  and  68®  V;.  (Fig.  2  a  and  b) 
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2.  Barotropic  .viodels 

The  first  attempt  to  explain  meanders  of  ocean  currents  from 
hydrodynamic  theories  was  almost  simultaneously  made  by  Haurwitz  and 
Panofsky  (1950)  for  the  Gulf  Stream  and  by  Ichiye  (195Q  for  the  Kuroshio. 

Both  theories  are  based  on  the  model  of  waves  along  the  boundaries  of  a 
jet  on  a  barotropic  rotating  fluid,  either  with  a  discontinuous  horizontal 
shear  (Haurwitz  and  Panofsky)  or  with  a  discontinuous  horizontal  velocity 
(Ichiye).  Since  the  velocity  profile  assumed  by  the  former  authors  is  more 
realistic,  their  result  is  discussed  here  and  referred  to  as  the  H-P  theory. 

The  basis  for  the  H-P  theory  is  observed  facts  that  the  ampli- 
tude  of  Gulf  Stream  meanders  begins  to  increase  after  the  Stream  leaves 
the  shelf  north  of  Cape  Hatteras,  and  that  on  the  coastal  side  of  the  Stream 
the  velocity  increases  to  its  maximum  considerably  stronger  than  decreases 
on  the  other  side,  forming  the  "shear  zone”  on  the  former  and  the  "transition 
zone"  on  the  latter  side.  The  theory  was  developed  in  order  to  explain  that 
the  effect  of  the  coastal  barrier  is  to  stabilize  lateral  waves  which  are  per¬ 
turbations  of  the  shear  zone  on  the  edge  of  the  stream.  In  the  H-P  theory, 
the  motion  is  assumed  to  be  independent  of  depth  because  the  shear  zone  is 
so  steep  that  the  perturbations  are  essentially  horizontal  and  thus,  the 
stabilizing  effect  of  gravity  may  be  neglected.  The  effect  of  the  earth's 
rotation  is  considered  but  not  its  variation  with  latitude  since  the  curvatures 
of  the  velocity  profile  across  the  current  are  large  compared  with  the  latter. 

The  equations  of  motion  and  continuity  for  the  perturbations  super¬ 
posed  on  the  horizontal  current  U  in  the  x-  direction  become  (after  omission 
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of  second  order  terms) 


9u/  9t  +  U8u/3x=  (f  -  9  U/  8y)  -  {  1/p  )  9p/9x 


(1) 


8v/9t+  U9v/  9x  =  -  fu  -  (  1/  p  )  9  p/  9  y 


(2) 


9u/  9  x  +  9  v  /  9y  =  0 


(3) 


It  is  assumed  that  the  basic  current  U  is  geostrophic  and  thus,  variable  in  a 
y-  direction  only. 

The  boundary  conditions  are  that  v  =  0  at  rigid  boundaries  and 
that  at  internal  boundaries  where  U  is  discontinuous 


(9  /  9t  +  U9/9y)  (p  -  p'  )  +  v9  (P  -  P'  )/  9y  =  0 


(4) 


and  an  analogous  equation  in  U'  and  v1  where  a  prime  indicates  quantities 
on  the  other  side  of  the  boundary  and  P  is  the  pressure  for  the  basic  current. 
(When  U  is  continuous,  v  and  p  must  also  be  continuous  across  the  internal 
boundary).  All  the  perturbation  quantities  are  assumed  to  be  of  form 


(u,  v,  p/  p  )  =  ( A,  B,  C  )  exp  Z-rr(±y+x-ct)/  L 


where  L  is  the  wave  length  and  c  the  wave  velocity.  When  the  wave  velocity 
is  complex,  the  real  part  gives  the  speed  of  propagation  and  the  imaginary  part 
c^  gives  the  rate  of  decrease  or  increase  of  the  amplitude  according  to  its 
negative  or  positive  sign.  In  different  zones  with  a  uniform  and  different  basic 
current  U,  the  constants  A,  B  and  C  can  be  determined  as  a  function  of 
one  of  them  from  equations  of  motion  (1)  and  (Z)  and  boundary  conditions  and 
the  last  boundary  condition  may  be  satisfied  provided  that  a  relation  (frequency 
equation)  exists  between  wave  speed  and  wave  length.  For  examole,  if  there 
are  two  zones  of  different  basic  currents,  the  number  of  constants  correspond¬ 
ing  to  A,  B  and  C  is  six  and  there  are  four  equations  of  motion  and  three 


~9~ 

boundary  conditions. 

Three  type*  9 i  velocity  profile*  for  the  basic,  current  which  is 

.  1 

K'  .  *  \  V 

bounded  by  solid  bar riera  on  both  sides  were  considered  in  the  H-P  theory. 

i  A*  J  *  .  •*»  't  s 

V?  *  ft,  *  * 

Model  1:  'f'he  sfmplesf  rnodef  consists  of  two  bounded  zones  with  speeds  of 

\  V  .  N  .'iv  ’ 

the  basic  current  aud  U2  ^ospqcfively,  of  width  d^  and  d2  (Fig,  3a). 

The  wave  velocity  is  giveh  by 

cslUjaj  +  U2a2)/(a^+a2)  ±  i  (  Uj  -  U2  )  /  N  a2/  aj  +  's/al/a2  )(5) 

Here  and  in  the  following  equations  an  coth  (2  d^/  L>.  It  is  seen  that  the 
waves  are  always  unstable.  Also,  the  imaginary  part  of  c  is  the  largest 
when  aj  and  a2  equal  unity.  This  condition  is  almost  satisfied  when  the 
discontinuity  occurs  more  than  half  a  wave  length  from  both  boundaries.  When 
dj  or  d^  is  much  smaller  than  L,  the  imaginary  part  tends  to  zero,  indi¬ 
cating  that  the  stability  increases  when  the  velocity  discontinuity  approaches 
either  boundary. 


Model  2:  For  the  model  which  consists  of  the  two  currents  with  speeds  Uj 
and  separated  by  a  shear  zone  with  the  current  varying  linearly  with  y 

(Fig.  3-b),  the  wave  speed  is  given  by 

c  =  (U1  +  U3)/2  +  (U3-U1)L(N-R)/4Trd2M  (6) 


+  (  U3  -  U ,  )  2. 


2 


1  -  (N  +  R)  L /  M  it  d,  +!L2  (N-R)2/  M  +  4  ^ /4ir2d?2  M* 


J 


where 


M  =  1  +  a  j  a^  +  aj  a3  +  a2  a3,  N  =  a^  +  a3  and  Rsaj+a^. 

For  long  waves,  for  which  a2  »  aq  and  a3  »  a^,  the  expression  under 

C  2  2,  Z 

the  square  root  becomes  11  -  J_  tt  +  (  d3  /  d2  )/  it  L  -  d3  LI  (tt  d2  )  which  is 

L  2  J 


1 


-  ^ 


••10'’ 


negative  in  most  cases,  excepting  if  the  shear  zone  is  so  close  to  one 
solid  boundary  tha.t  d^/  d2  is  nearly  equal  to  zero.  This  again  indicates 
that  the  solid  boundary  stabilizes  the  perturbations  on  the  shear  zone. 


Model  3:  A  better  approximation  to  the  jetlike  velocity  profile  of  the  Gulf 
Stream  by  four  zones  where  the  current  speeds  in  the  two  outer  zones  are 
constant  and  between  them  lie  the  shear  and  transition  zone  in  which  the 


velocity  gradients  are  constants,  as  shown  in  Fig.  3-c.  The  frequency 
equation  is  a  cubic  of  velocity  c  but  the  coefficients  are  quite  complicated 
However,  in  the  case  that  Uj  =  U^j  a  condition  which  seems  to  be  reason¬ 
able  since  the  perturbations  to  be  studied  are  the  meanders  along  the  Gulf 
Stream,  the  frequency  equation  becomes  much  simpler  and  is  given  by 


c^A-cg2  A  +  S3(B-E)  +  S2(D-E)!  +Cj 


/ 


A 


S3  E  +  S2  D  -  S3  J 


-  S2  I-S2S3(  J  +  I-H) 


+  S2S3(S2  +  S3-H)  =  0 


where 


cL  =  (  c  -  Uj  )  /(U2  3  -  Ux) 

A  =  (  aj  +a2  )  (  1  +  a3a4)  +  (  a3  +  a4 )  (  1  +  aj  a^  ) 

B  =  1  +  ax  a2  +  a,  a3  +  a2a3 

D  =  1  +  a0  a0  +  a0  a  .  +  a0  a  . 

2  3  2  4  3  4 

£  =  (aL  +a2)  (a3+  a4) 

J  =  aj  +  a2,  H  =  a2  +  a3,  I  =  a3  +  a4 
S2  =  L  /  2  tt  d2  and  S3  =  L  /  2  ir  d3 


(7) 
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Equation  (5)  indicates  that  depends  not  on  the  numerical  values  of  the 
basic  current,  but  on  the  geometry  of  the  zones  only.  Since  d^  may  be 
considered  as  the  distance  between  the  Gulf  Stream  and  European  conti¬ 
nent,  it  is  safely  assumed  that  d^  »  L  or  a^  =  1.  The  unstable  waves 
can  develop  only  if  the  equation  has  conjugate  complex  roots  which  are 
possible  when  its  discriminant  is  positive.  The  values  of  the  discriminant 
for  different  values  of  d^  /  L  are  given  in  the  following  table  in  a  model 
with  d^/ L  =  1/44  and  d^/L  =  1  / 11*4  corresponding  roughly  to  the  Gulf 
Stream: 

d4/L  1/63  1/630  1/4400  1/6300  0 

Disc.  1.523  0.244  0.024  0.006  -0.0005  -0.0030 

This  table  suggests  that  the  decrease  of  the  distance  between  the  shear  zone 
and  the  continent  increases  the  stability  of  the  current.  A  more  instructive 
presentation  is  given  by  Fig.  4  v/hich  shows  curves  of  conjugate  roots  versus 
wave  length  L  of  Equation  (5)  for  different  values  of  d^  (Ichiye,  1955,  1956). 
In  this  model,  d^  =  100,  d^  =  20  and  d^  =  20  (curve  I)  and  200  (curve  II), 
all  in  nautical  miles.  When  the  wave  length  is  below  a  critical  value,  the 
conjugate  roots  of  (5)  have  two  real  values  and  the  waves  are  stable.  When 
the  wave  length  exceeds  the  critical  one,  the  roots  become  complex  and  the 
real  and  imaginary  part  are  plotted  as  separate  curves  in  this  figure.  The 
real  part  gives  the  propagation  speed  of  the  disturbance  while  the  imaginary 
part  represents  the  rate  of  growth  of  the  amplitude.  This  figure  indicates 
that  the  longer  the  distance  between  the  coast  and  the  shear  zone,  the  smaller 
the  critical  wave  length  at  which  the  instability  occurs  and  the  greater  the 


{ 1 
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growth  rate  of  the  amplitude  of  unstable  waves.  For  comparison,  by  use 
of  a  theory  developed  by  Kuo  (1949),  the  real  and  imaginary  parts  of  the 
dimensionless  wave  velocity  (c  -  UQ/  -  UQ)  were  computed  with  a 
model  of  planetary  waves  in  a  zonal  jet  which  has  a  horizontal  velocity 
profile,  as  shown  in  Fig.  3-d  (Ichiye,  1956).  In  this  calculation,  the  width 
of  the  jet  was  taken  as  120  n.  miles  and  the  velocities  Um  and  UQ  were 
taken  1.  1  and  -0.1  knots,  respectively.  It  is  seen  that  the  critical  wave 
length  for  this  model  (planetary  waves)  is  much  larger  than  for  the  H-P 
model  (boundary  waves)  and  also,  the  imaginary  part  of  the  wave  velocity 
or  growth  rate  of  the  amplitude  of  unstable  waves  is  much  smaller.  However, 
the  growth  velocity  for  the  planetary  wave  model  has  maximum  value  at  the 
wave  length  of  about  250  n.  miles,  corresponding  to  the  prevalent  meanders 
in  the  Gulf  Stream  and  the  Kuroshio,  while  the  one  for  the  H-P  model  increases 


with  the  wave  length  and  thus,  it  has  no  prevalent  wave  length. 
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3.  Baroclinic  Waves  (case  of  no  horizontal  shear) 

Stommel  (1953)  considered  the  long  waves  in  a  stratified  ocean 
in  which  horizontal  divergence  was  dynamically  important.  In  order  to 
simplify  the  analysis,  he  abandoned  the  horizontal  shear  of  the  basic  cur¬ 
rent  and  assumed  a  two-layer  ocean,  in  which  the  lower  layer  is  suffi¬ 
ciently  deep  so  that  there  are  no  horizontal  pressure  gradient  and  current 
there.  The  geostrophic  current  U  in  the  upper  layer  flows  in  the  x- 
direction.  A  cross-stream  pressure  gradient  due  to  this  current  is  given 
by 

fU  =  -g3H/9y  =  -  g  (3D/3y)  Up/p  ) 
in  which  H  is  the  elevation  of  the  free  surface,  D  is  the  thickness  of  the 
upper  layer  and  £p  is  the  difference  in  density  of  the  upper  and  lower 
layer. 

When  small  perturbations  of  u,  v  in  velocity  and  h  in  H  are 
assumed  to  be  independent  of  y,  the  perturbation  equations  of  motion  and 
continuity  may  be  written  in  the  form: 

du/dt-fv  =  -g9h/9x;  dv/dt+fu  =  0  (1)(2) 

and 

(  1  +  p  /  £  p  )  dh  /  dt  +  D9u/9x  +  v9D/9y  =  0  (3) 

where  d/dt  is  the  operator  (9  /  9t  +  U9/9x).  If  the  perturbations  are  all 
in  the  form  expi  (  kx  -  <rt),  the  frequency  equation  becomes 

k2U2  (1  -  b)3  -  (f2  +  g'k2D)  (1-b)  +  f2  =  0  (4) 

where  g'  =  g  l  p/p»  b  =  c/U  =  <r/kU  =  U+  b"i  and  b'  and  b"  are  real. 
Equation  (4)  can  be  rewritten  as 

Y3  +  2  -  6  Y  =  0 


where 


(5) 
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6  =  (2  f  /k  U) 2/ 3  (1  +  g'k2  D  fZ) 

Y  =  Z  (2  f/kU)"2^3  (1  «  b) 

The  roots  of  Equation  (5)  are  all  real  if  6  >  3  and  have  conjugate  complex 

values  if  6  <  3.  As  in  the  previous  case,  real  and  conjugate  complex  roots 

correspond  to  stable  and  unstable  waves,  respectively.  The  stability 

criterion  6  >  3  becomes  UZ  <  g'D.  When  UZ  =  g'D,  the  instability 

begins  with  vanishing  b'  or  stationary  waves,  wave  number  of  which  is 

k  =  f  / n/  2  U.  Unstable  waves  have  a  narrow  range  of  wave  numbers  about 

this  value  for  slightly  larger  U  than  the  critical  value.  'When  numerical 

values  of  D  =  200  m,  U  =  2  m  /sec  and  Ap/p  =  2  x  10"  are  used,  the 

critical  wave  length  becomes  180  km.,  which  corresponds  to  that  of  the  large 

stationary  meander  observed  in  Operation  Cabot.  Rossby  (1951)  also  showed 

1/2 

that  the  velocity  of  the  Gulf  Stream  approaches  a  critical  value  (g'D)  in 
a  steady  state  due  to  vertical  concentration  of  momentum,  and  that  the  steady 
state  may  be  broken  down  if  the  velocity  exceeds  this  value.  However,  there 
is  no  indication  that  the  depth  D  of  the  Gulf  Stream  decreases  sharply  be¬ 
tween  Cape  Hatteras  and  the  tail  of  the  Grand  Banks  and  thus,  this  theory 
does  not  explain  the  sudden  increase  of  the  meanders  east  of  about  65°  W, 
as  described  in  the  previous  section, 

Iida  and  Ichiye  (1963)  showed  that  there  is  no  unstable  regime  for 
the  perturbations  in  the  uniform  geostrophic  current  of  a  two-layer  ocean 
when  the  linear  increase  of  D  with  the  lateral  distance  y  is  taken  into 
account.  Stommel  assumed  that  the  thickness  of  the  upper  layer  D  is  con¬ 
stant  and  also  its  lateral  gradient  9D/3y  is  finite  and  constant.  This 


-15- 


assumption  may  be  permissible  in  the  central  part  of  the  ocean,  but  incom¬ 
patible  along  the  edge  of  the  Gulf  Stream,  where  the  pronounced  meanders 
were  observed.  When  there  is  a  uniform  current  U  in  the  upper  layer,  the 
depth  D  decreases  linearly  with  y  to  the  left  of  the  current  (facing  down¬ 
stream)  and  thus,  the  interface  between  the  two  layers  intersects  the  surface 
(Fig.  5).  When  this  intersection  is  taken  as  the  x-axis,  the  two  layer  ocean 
extends  to  the  negative  y-plane  and  D  is  expressed  by  -  sy,  where 
-  s  =  f  U/  g'  ,  Since  the  perturbations  u,  v,  and  h  are  dependent  on  y  as 
well  as  x,  the  righthand  side  of  Equation  (Z)  is  replaced  by  -  g  3  h/  3y.  The 
elimination  of  u  and  v  in  Equations  (1)  to  (3)  leads  to  the  equation  for  h: 


d  (  D  dh/  dy )  /  dy  +  (fks/ffj  -  k^D  -  /  g' s  )  h  =  0  (6) 

2  2 

where  =  cr  -  kU  and  A  =  f  -  .  Stommel's  frequency  equation 

may  be  derived  by  equating  the  coefficient  of  h  in  Equation  (6)  to  zero.  When 


D  -  sty  is  substituted  in  this  equation,  we  have 


d  (y  dh  /d y)  - 


k/u>  +  f  (1  -  to2)  /  U  +  k2y 


h  =  0 


(?) 


where  to  =  <r  ^  /  f. 

In  order  that  the  normal  mode  of  h  which  is  finite  at  y  =  0  and 


c.*~,  may  exist,  to  should  satisfy 

k  to"1  +  f  (  1-co2)  U-l  =  -  (2n  +  1)  |fc\ 


(8) 


in  which  n  is  positive  integers.  Equation  (8)  can  be  rewritten  as 

W3  +  2  =  6'  W 

t 

where  6'  is  given  by 

6'  =  (  2/ k')~2/3  (2n  +  1  )  k'  +  1 

and  k'  »  kU/f  and  W  =  -  to  (k*  /  2)  ~  ^  3 .  As  before,  Equation  (9) 


(9) 


(10) 
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has  conjugate  complex  roots  when  6'  <  3.  However,  it  is  seen  that  for 

all  positive  k'f  6' >.3,  where  the  equality  occurs  for  k'  =  2  and  n  =  0. 
It  is  also  proven  that  the  frequency  equation  has  only  real  roots  when 
there  is  a  constant  geostrophic  current  U'  in  the  lower  layer.  In  this 
case,  k'  in  Equation  (10)  should  be  replaced  with  (1  -  U '/  U )  k' ,  Since 
the  latter  quantity  is  always  positive  from  geostrophic  assumption,  the 
same  discussion  as  in  the  case  of  motionless  lower  layer  is  valid.  It  is 
interesting  to  note  that  Equation  (9)  becomes  equal  to  the  frequency  equa¬ 
tion  (30)  in  Reid's  paper  (1958),  when  u>  and  k'  are  replaced  by  -co/f 
and  g'sk/f,  respectively.  Therefore,  mathematical  characteristics  of 
roots  of  Equation  (9)  are  the  same  as  those  for  the  barotropic  waves 
propagating  along  a  coast  in  a  sea  with  a  depth  linearly  increasing  from 
the  coast. 

Although  Stommel's  original  solution  has  an  interesting 
feature  for  stability  of  a  baroclinic  current,  it  has  some  inconsistency. 

In  order  to  discuss  the  stability  problem,  we  must  take  a  lateral  shear 
of  the  current  into  consideration. 
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4„  Baroclinic  Waves  (model  with  horizontal  shear) 

Iwata  (1961)  and  Duxbury  (1963)  treated  the  perturbations  in 
a  two -layer  ocean  in  which  the  geostrophic  current  in  the  upper  layer  is 
not  constant,  as  in  Stommel's  case,  but  is  variable  along  the  lateral 
direction,  Iv/ata  discussed  the  solutions  for  the  cases,  in  which  a  horizontal 
velocity  distribution  is  given  by  different  profiles,  including  a  simple  jet 
of  a  uniform  velocity  in  an  infinite  motionless  fluid  and  a  uniform  mean 
current  in  a  half  plane  with  another  motionless  half  plane.  He  employed 
the  assumption  of  the  geostrophic  perturbation  velocities.  This  assump¬ 
tion  is  valid  when  jcr  -  kU|  <<  f,  where  cr  and  k  are  frequency  and 
wave  number  of  the  perturbations,  respectively.  On  the  other  hand, 

Duxbury  treated  only  the  case  of  the  uniform  mean  flow  in  a  half  plane 
bounded  by  the  motionless  half  plane  but  without  the  assumption  of  geo¬ 
strophic  perturbation  velocities,  thus  covering  more  general  range  of 
the  frequency  than  Iwa„ta. 

The  equations  of  motion  and  continuity  for  the  perturbations  in 
a  two-layer  ocean  with  the  motionless  lower  layer  of  an  infinite  depth  are 
almost  the  same  as  Equations  (1)  to  (3)  in  Section  3,  except  -g  3  h/ 8  y 
replacing  zero  in  the  righthand  side  of  (2).  (Hereafter  an  equation  in 
other  sections  will  be  distinguished  by  the  Roman  numeral  denoting  the 
section;  for  example,  III-l).  However,  in  most  of  the  recent  literature, 
the  thickness  of  the  upper  layer  instead  of  the  elevation  of  free  surface  is 
used  as  a  variable  and  thus,  the  perturbation  equations  are  derived  here 
again. 


■ 


-  l  a. 


If  V  (x,  y,  t)  denotes  the  total  horizontal  components  of 
velocity  in  a  two-layer  ocean  with  the  still  lower  layer  of  an  infinite  depth 
and  H  (x,  y,  t)  the  thickness  of  the  upper  layer,  the  momentum  and  continuity 
equations  in  the  vectorial  form  become 

av/9t  +  V  •  vv  +  zi  X  $  =  -  g»  VH  (1) 

8H/a  t  +  v  V  H  =  0  (2) 

in  which  g1  =  g  l.  p  /p  is  the  "reduced"  gravity  and  oJ  is  the  angular 
velocity  of  the  earth's  rotation.  The  velocity  can  be  resolved  into  a  mean 
current  U(y)  directed  along  the  x-  axis  and  perturbations  whose  components 
are  u  and  v.  With  a  similar  decomposition  of  H  =  D  +  h,  the  linearized 
perturbation  equations  become 

du/dt+(dU/dy-f)v+g'8h/3x  =  0  (3) 

dv/dt+fu+g'8h/8y  =  0  (4) 

dh/dt  +  8  (Du)/ax  +  8(Dv)/8y  =  0  (5) 

where  (d/dt)  indicates  8/t+U8/8x.  Notice  that  here  h  is  not  the 
surface  elevation  but  the  perturbation  of  the  thickness  of  the  upper  layer. 

The  boundary  condition  along  the  lateral  boundary  between  discontinuous 
mean  flows  corresponding  to  (II-4)  becomes,  with  neglect  of  second  order 
terms, 

(8/8t  +  Uj8  /  8x)  'hi  -  h2  )  +  v  d  (Dj  -  D2)  /  dy  =  0  (6) 

for  j  =  1  and  2,  where  the  subscripts  denote  regions  with  different  mean 
flows.  As  a  first  approximation,  Equation  (6)  should  be  satisfied  at 
y  =  yQ  (constant),  which  is  the  y  coordinate  of  the  undisturbed  boundary 
Df  the  two  mean  flows  and  U^.  The  lateral  excursion  of  the  boundary 

*  can  be  expressed  in  terms  of  hj  and  h^  from  the  condition  that 


X 

. 
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Hj  =  at  y  =  y0  +  £  .  Thus,  the  first  approximation  for  £,  is  given  by 

4  =  (hj  -  h2)  jd  (D2  -  Dl  )  /dyj  _1 

where  h,  and  h_  are  evaluated  at  y  =  y  . 

12  o 

In  Duxbury's  model  which  is  shown  in  Fig,  6a,  the  uniform 
current  U  extends  in  the  region  1  (y  <  0),  while  the  region  2  (y  >0) 
is  motionless.  The  thickness  D  in  the  region  1  and  2  is  given  by 


(7) 


DQ  -  sy  and  D0,  respectively, 
of  the  form 


The  variables  u.,  v.  and  h.  are 

J  J  J 


(uj*  vj-  V = 


X.  (y)>  Y.  (y)  ,  z  (y)  ;  exp  i  (kx  -  cr  t) 
J  J  J  i 


Elimination  of  Xj  and  Yj  in  Equations  (3)  to  (5)  leads  to 


(8) 


(D  -  sy)  d2Zj/  dy2  +  s  d  Zj/dy  +  (K  -  k2DQ  +  k2  sy)  Zj  =  0  (9) 


where 


d2  Z2/dy2  +  mZZ2  =  0 


K  =  (7  -  f2)  IgT1  - 


f  s  k 


-1 


u> 


m  = 


(>2  -  f2  )  (g’  D0>-‘  -  k2 


1/2 


U) 


=  0-  -  kU 


(10) 

(11) 

(12) 

(13) 


If  the  perturbation  velocities  are  assumed  to  be  geostrophic,  co  and  <r 

are  neglected  against  f2  in  Equations  (11)  and  (12),  The  boundary  condition 

(6)  which  couples  functions  Zj  and  Z^  becomes 

iu.(Z,  -ZJ  +  sY.  =0 
J  1  L  j 


(14) 


wnere 


co  i  —  oj  ,  u>t“  o"  and 


Y.  =  -  g'i  (fkZj  +  oojdZj/dy)  (/.  -  fA) 


2.-1 


(15) 


The  solution  of  (9)  which  satisfies  the  condition  of  finiteness  at 


y  =  -  Oo  is  given  by 


7  =  Bj  V»  (  a,  z) 


(16) 


No¬ 

where 

W  (  a  ,  z)  =  exp  (~  z/2)  ^ (a  ,  1;  z)  (1 

in  which 

a  =  (  1/2)  (  1  -  K/  sk)  (If 

z  =  2k  (CQ  -  sy)  /  s  (19a 

and  St  (a,  1;  z)  is  the  Tricomi  function  (Erdelyi  et  al,  1953).  Equation 
(16)  holds  for  z  >  z^, ,  where 

zo  =  2kDo/s  <19b' 

For  y  >  0,  the  solution  of  (10)  which  satisfies  the  condition 

of  finiteness  at  y  =  «>„■  is  given  by 

exp  (~  y)  (20) 

2 

when  m  is  negative  (=  -  vl),  or  when 

a1  >  f2  +  g’ L0k2  (21) 

In  this  case,  the  perturbations  in  region  1  become  the  Kelvin-type  boundary 

waves.  When  is  positive,  instead  of  (20),  we  have 

Z£  =  A'^  cos  my  +  B^'  sin  in  y  (22) 

2 

Duxbury  did  not  discuss  the  case  of  positive  m  .  However,  in  this  case 
when  there  is  a  coast  at  a  finite  distance  from  y  =  0,  can  be  expressed 

with  from  the  boundary  condition  that  v  =  0  at  the  coast.  Then,  the 

subsequent  process  for  obtaining  the  frequency  equation  is  the  same  as  for 
solution  (20).  If  there  is  no  coast  on  the  side  of  positive  y,  there  is  no 
relationship  between  a  and  k,  except  the  inequality  (21)  with  the  reverse 
sign.  Then,  the  problem  is  to  determine  the  relationships  among  the  inci¬ 
dent  (oblique)  waves  and  the  reflected  waves  in  region  1  and  the  refracted 
waves  in  region  2  (Ichive,  1963). 

Solutions  (16)  and  (20)  should  satisfy  two  boundary  conditions 


.  •>.  •  '  1 
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(14)  at  y=  0.  Then,  elimination  of  Bj  and  leads  to  the  frequency- 

equation  ^ 

w  (  co2  -  f2  )  -  g'  sk  )f  -  2W  '  (a,  zq)  /  Vr  (a  ,  zq)  | 


cr  ( cr2™f2)  +  g*  sk  (f  -  <r  lA.  /  k ) 


=  (  oo2  -  f2)(<r2  -  f2)  (23, 


where  W  (a,  z)  denotes  dV// dz.  Introducing  the  dimensionless  vari¬ 


ables 

where 


P  =  w/f,  j..  =  Uk/f  =g'skf"2  =  F2  zq/2 
F  =  U  /  (g‘  D0)~ 1/2 


(24  a,  t 
(24c) 


Equation  (23)  can  be  rendered  in  the  form 


^  ^  (CE  -  1)  +  (E  -  v)  (1/+  (  |/+  |,)2  -  1( 


in  which 


I  l 

4 


b  s 


-  <v+  L)2  -  l  ( 


«■  =(l/2) 


1  /  2  1\  **  1  —1 
1  -  (  ^  -  1  )  p  +  V 


C  —  2  cl  —  1 


E(  a,  zQ)  =  V:  (a,  zq)  /  2  (a  ,  zQ) 


(25) 

(26) 

(27) 

(28) 


The  parameter  y  is  a  relative  "Doppler"  frequency,  p  is  a  relative 
longitudinal  wave  number  and  F  is  the  internal  Froude  number.  Equation 
(25)  is  an  eight-degree  polynomial  in  p,  in  which  the  coefficients  depend 
upon  a  ,  p  and  F,  and  (26)  is  a  third-degree  polynomial  in  p.  A  practi¬ 
cal  method  to  obtain  the  relation  between  p  and  p  is  first  to  determine 
graphs  of  p  versus  p  for  selected  values  of  a  from  Equation  (26), 

Then,  from  Equation  (25)  graphs  of  p  versus  a  are  determined  for 

? 

particular  values  of  F  and  a,  using  the  computer.  Each  graph  was 


-2  2- 


divided  along  its  p.  -  axis  into  zones  of  different  mode  number  (n),  which 
was  determined  by  the  number  of  zeros  in  argument  z  of  W  (  a  ,  z)  or  a 
transverse  wave  number.  Then,  the  graphs  obtained  from  Equations  (25) 
and  (26)  are  superposed  to  determine  the  admissible  values  of  v  and  a 
for  common  p  and  F,  corresponding  to  each  mode  number  n.  Further 
V  and  p,  must  satisfy  condition  (21)  or  its  dimensionless  form 


(  P  +  p)2  -  1  <  (p  /  F) 


(29) 


One  of  the  jyj-  p  graphs  for  F^  =  10  is  shown  in  Fig.  6b.  These  graphs 
indicate  that  there  are  two  classes  of  waves:  inertio  -gravitational  waves 
with  \i-  >  1  (orjoji  >  f),  and  quasi-geostrophic  ones  with  |  i'  <  1  (or 
|  u>  |  <  f)»  each  class  having  different  mode  numbers.  Only  the  zero-order 

mode  behaves  like  the  quasi-geostrophic  class  for  small  p,.  The  inertio- 
gravitational  waves  travel  upstream  while  the  quasi-geostrophic  waves 
travel  downstream.  The  former  has  small  wave  lengths  and  it  appears  to  be 
irrelevant  to  the  meanders.  The  physical  parameters  for  the  Gulf  Stream, 
such  as,  Dq  =  100  m,  s=  1.5x10"*,  6  p  /  p  =1,1x10"*  ,  U  =  1.8  m/sec 

leads  to  F  =  3.  Duxbury's  analysis  indicates  that  there  are  stable  waves 
of  the  quasi-geostrophic  class  of  wave  lengths  and  speeds  comparable  to 
those  of  the  Gulf  Stream  meanders  when  the  internal  Froude  number  (F) 
exceeds  unity,  unlike  the  result  of  Stommel  (1963).  Further,  Duxbury 
considered  that  the  Gulf  Stream  meanders  could  be  generated  by  the  stable 


modes  from  three  reasons:  (1)  The  disturbances  composed  of  many  modes 
cf  dispersive  character  would  grow  in  amplitude  downstream  by  manner  of 


-2  3- 


phase  reinforcement  (2)  An  external  energ-^Like  wind  torque  would  cause 
resonance  coupling  with  the  stable  waves  which  are  actually  of  the  normal 
modes  of  the  system  and  (3)  Variations  of  Dq  and  U  may  produce 
amplification  analogous  to  the  amplification  of  swell  propagating  into  shal¬ 
low  water.  However,  his  analysis  does  not  prove  non-existence  of  the 
unstable  waves  (or  complex  roots  of  the  frequency  equation)  in  the  range 
of  wave  lengths  and  propagation  velocity  comparable  to  the  actual  meanders. 
The  first  two  possibilities  of  exciting  meanders  by  stable  modes  seem  to  be 
unreasonable  since  occurrence  of  a  localized  disturbance  and  an  external 
energy  source  is  random  and  yet  the  actual  meanders  are  rather  regular  in 
their  wave  characteristics  and  their  geographical  locations.  Also,  without 
evaluation  of  the  solutions  corresponding  to  the  initial  condition  of  the 

localized  disturbance,  the  external  force  or  the  variable  D  and  U, 

o 

the  three  possibilities  described  by  Duxbury  are  mere  speculations.  As 
far  as  the  mechanism  of  the  Gulf  Stream  meanders  is  concerned,  Duxbury' s 
analysis  does  not  give  any  light  on  the  problem. 

Iwata  (1961)  derived  the  frequency  equation  si^Dar  to,  but 
simpler  than  (23)  under  the  assumption  of  geostrophic  perturbations  for  the 
same  profile  of  the  mean  flow  as  above.  He  did  not  solve  the  frequency 
equation  numerically  nor  treated  it  analytically,  in  general  termr.  However, 


he  could  show  that  when  zQ  »  1,  the  frequency  equation  becomes 


2  .  ^2  ,,  .  -1  -1,,-!.  -1 


k  -  j  p 


+  F£  (  1  +  |*  |*1  )/2|i  +  n  (.J  =0 

J 

from  the  asymptotic  expansion  of  V?  (  a  ,  z)  for  large  z,  where 
L  =  cr/(kU-cr)  =  -(r/yf 


(30) 


(31) 


M  =  b  +  F' 


(32) 
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Also,  when  zQ  <  <  1,  the  frequency  equation  becomes 

2i2  +  FZL  -  F  =  0  |33) 

It  is  easily  seen  that  Equation  (30)  has  conjugate  complex  roots,  while 

I 

(33)  has  real  roots  only.  The  definition  of  zQ  given  by  (19b)  yields 

zq  =  1.3  x  10  k  (k  in  m  )  from  the  physical  constants  corresponding  to 

the  Gulf  Stream  (D  =  100  m,  s  =  1,5  x  10"^  m ** ^ ) .  Therefore,  fc  =1 

°  o 

yields  the  perturbation  wave  length  of  84  km.  This  suggests  that  in  the 
range  of  wave  lengths  comparable  to  the  real  meanders,  the  condition 
zQ  <<  1  is  satisfied  and  thus,  the  perturbation  is  stable,  as  discussed 

■  I !  1 

by  Iida  and  Ichiye  (1963),  In  another  example,  Iwata  (1961)  treated  the 

(it  f 

case  of  the  jetlike  mean  flow  with  a  constant  velocity  and  a  finite  width, 
as  shown  in  Fig.  6c.  V  hen  the  thickness  of  the  upper  layer  of  the  jet  is 

l  i  kI  * 

a  linear  function  of  y,  satisfying  the  geostrophic  relation  for  the  mean 

■ » ’  : 

t  W 

flow,  the  solution  for  the  jet  region  corresponding  to  Equation  (16)  includes 

i  *.  ■  Jli  ht 

the  Kummer  function  £  (a,  c;  z)  (Erdelyi  et  al,  1953)  as  well  as  the 

t  w  -i 

Tricomi  function  and  the  analytical  treatment  of  the  resultant  frequency 

4 

H 

equation  becomes  intangible.  Iwata  alternatively  assumed  that  the  depth 
of  the  upper  layer  of  the  jet  is  constant  on  the  basis  of  smallness  of  its 
velocity  and  obtained  the  frequency  equation  which  gives  unstable  waves, 

I 

>  il 

under  the  assumptions  that  the  wave  lengths  of  the  disturbances  are  much 

1  ;  f  !' 

smaller  than  300  times  the  width  of  the  jet.  His  frequency  equation  is 
only  an  approximation,  but  it  is  certain  that  there  is  a  range  of  unstable 

l  4'\ 

waves  with  wave  lengths  comparable  to  the  real  meanders.  Also,  validity 
of  his  assumption  of  the  uniform  depth  of  the  jet  must  be  checked  by  scaling 
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basic  equations  using  characteristic  quantities  as  will  be  discussed  below. 

However,  for  wave  lengths  much  larger  than  the  width  of  the  jet,  this 

assumption  might  be  valid.  This  is  similar  to  a  case  discussed  by  Ichiye 

(1951)  on  propagation  of  long  waves  in  a  channel  of  variable  depth  and  width 

connected  with  another  uniform  channel.  It  was  shown  that  the  variable 

channel  can  be  replaced  with  a  uniform  channel  with  averaged  values  of 

depth  and  width  of  the  prototype  channel  when  the  wave  length  is  sufficiently 

long.  It  will  be  worthwhile  to  test  the  applicability  of  this  result  to  the 

quasi-geostrophic  waves  in  a  geostrophic  jet  considered  here, 

Yasui  (1961)  treated  the  perturbations  in  the  same  two-layer 

ocean  as  in  Iwata's  and  Duxbury's  model.  He  obtained  an  exact  solution  of 

the  equation  of  the  surface  elevation  in  terms  of  a  hyper-geometric  function 

2 

for  the  distribution  of  the  mean  flow  given  by  U  =  UQ  sech  (y/L)  (Bickley 
jet),  under  the  assumptions  that 

<r  »  k  U  and  \cr2-f2]»  f)dU/dy|  (34  a,  b) 

The  frequency  equation  in  this  model  is  given  by  _1 

d  i  (^2  +  F2  <  1  -  K  2  ) 

’ 

where  n  is  non-negative  integers.  The  dimensionless  parameters  in  this 
equation  as  defined  by 

jjL  =  UQ k  /  f;  k  =  o- /  f  (3  6  a, 

d  =  2Lf/UQl  d*  =  LfDj/U  DQ  (37  a, 1 

while  the  thickness  of  the  upper  layer  is  given  by 

D  =  D0  +  Dx  tan  h  (y/L)  (38) 

Equation  (35)  becomes  a  fifth  degree  polynomial  in  k  after  rationalization. 


1/2  =  1  -  2n  +  |  1  +  d‘(  2  +  p/  ?  ) 


1/2 


I 
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Yasui  did  not  discuss  the  equation  in  detail.  Also,  the  assumptions  (34  a,  b) 
may  exclude  the  waves  with  a  low  propagation  speed  corresponding  to  the 
meanders.  It  will  be  worthwhile  to  treat  analytically  the  case  of  the  Bickley 
jet  for  the  low  propagation  speed.  A  conclusion  immediately  obtained  from 
the  frequency  Equation  (35)  is  that  p,  k  should  be  positive  for  positive 
indicating  that  the  inertio- gravitational  waves  along  a  Eickley  jet  propagate 
in  an  opposite  direction  to  the  mean  flow,  as  in  the  case  of  a  shear  boundary 
discussed  by  Duxbury  (1963). 

Yasui  (1963)  also  treated  the  quasi-geostrophic  waves  in  a  jet 
shown  in  Fig.  6(c)  of  a  two-layer  ocean  by  assuming  the  depth  of  the  jet  to 
be  constant  in  a  way  similar  to  Iwata's  model  (1963).  He  applied  his  result 
to  fluctuations  of  subsurface  temperature  and  density  with  periods  longer 
than  a  week  observed  at  an  ocean  weather  station  T  at  29*  N  and  135®  E.  He 
speculated  that  the  fluctuations  of  periods  of  80  to  100  days  were  due  to  the 
stable  waves  of  a  quasi-geostrophic  mode  having  a  propagation  speed  close 


to  the  mean  flow. 


t  « 
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5.  Hydrodynamic  Stability  of  the  Gulf  Stream 

A  general  theory  of  hydrodynamic  stability  (Lin,  1955)  has 
been  applied  to  the  study  of  perturbations  in  the  atmospheric  jet  by  many 
authors  including  Kuo  (1949),  whose  method  of  studying  the  stability  of  a 
barctrcpic  flow  was  also  used  in  the  problem  of  ocean  current  meanders; 
for  instance  by  Ichiye  (1955),  as  discussed  before.  Stern  (1961)  and 
Pedlosky  (1963)  treated  the  stability  of  a  linear  flow  with  a  horizontal  shear 
in  a  two  layer  ocean. 

Stern  (1961)  discussed  the  energy  balance  of  a  baroclinic  jet 
in  a  two-layer  ocean.  Tee  momentum  and  continuity  equations  (III  -  1,2) 
are  valid,  if  the  hydrostatic  approximation  may  be  defended.  Then,  the 
two  layer  model  is  equivalent  to  a  single  layer  system  with  a  "reduced" 
gravity  and  with  the  depth  equaling  the  thickness  of  the  upper  layer.  -The 
condition  for  the  validity  of  the  hydrostatic  approximation  is  that  Dq/L 
<<  D^/L<<  1  when  R  £  1  and  F  4  1,  where  Dq,  and  L  are  the 

mean  depth  of  the  top  layer  and  bottom  layer  and  the  horizontal  scale  of 
perturbation,  respectively  and  R  and  F  are  the  stream  Rossby  number 
and  the  rotational  Froude  number,  respectively,  defined  by 

R  =  UQ/f  L  and  F  =  f2  L2/  g'  DQ  (  1  a,  b) 

Energy  equations  may  be  obtained  by  multiplying  the  vectorial 
form  combining  (III-3)  and  (III-4)  by  Jv  and  Equation  (III  —  5)  by  g'h, 
where  v  is  the  vector  of  the  perturbation  velocity  (u,  v)  and  U  and  D 
are  functions  of  y  only.  The  results  are 
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9E/3t+  U9E/8x+  Duv  (dU/dy)  +  g'Dv  •  Vh  =0  (2) 

9  C/3  t  +  U8Q/3x  +  g'h  V*  Dv  =  0  (3) 

where 

E  =  (  1/2)  D  Jvj  2  and  Q  =  (1/2)  g'h2  ( 4  a,  b) 

are  respectively  the  kinetic  and  potential  energy  of  the  perturbation  in  the 
undist  urbed  depth  D„  If  Equations  (2)  and  (3)  are  averr.ged  over  the  entire 
horizontal  plane,  the  second  terms  of  both  Equations  (2^  and  (3)  vanish, 
when  E  and  Q  are  finite,  because  U  is  the  function  of  y  only.  Also, 
the  average  of  the  last  term  of  Equation  (2)  equals,  but  has  a  reverse  sign 
to  the  average  of  the  last  term  of  (3),  bee?  :,se 


hV 


v  =  i  h  E  v 


n 


V 


Vh 


(5) 


from  Green's  theorem  where  is  the  velocity  component  normal  to  the 

boundar  y  and  the  brackets  indicate  the  integral  along  the  boundary  divided 

by  the  *otal  area.  The  quantity  in  the  brackets  vanishes,  when  the  v 

is  n 

\ 

vanishes  at  the  x~  boundari.es  ar.dxperiodic  in  x  far  downstreaxai.  The  right- 
hand  side  of  (5)  multiplied  by  g'  represents  the  work  done  against  the  gravi¬ 
tational  field.  Summing  averaged  Equations  (2)  and  {3},  we  have 


3(E  +  C)/3t=-Duv3U/3y  (6) 

which  indicates  that  the  rate  of  increase  of  total  energy  of  the  perturbation 
depends  on  the  Reynolds  stress  a.cting  on  the  horizontal  shear  of  the  mean 
flow.  Equation  (6)  implies  that  the  meander  theory  of  Stommel  (1953)  dis¬ 
cussed  in  Section  2  is  untenable,  for  in  his  model,  3  U/3  y  =  0  and  thus, 

E  +  Q  is  conserved.  Also,  this  equation  indicates  that  the  instability  is 
caused  only  through  the  shearing  mechanism  even  when  there  is  a  conversion 


. 
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of  mean  potential  energy  to  the  perturbation,  owing  to  the  hydrostatic 
assumption. 

Equations  (1IJ-1)  and  (III — Z )  imply  the  conservation  of  the 
vertical  component  of  potential  vorticity 

6(f+Z+4)H  /  6  t  =  0  (6) 


where  Z  and  are  the  relative  vorticity  of  the  mean  motion  and  the  per¬ 

turbation,  respectively.  H  is  the  total  depth  or  the  sum  of  mean  depth  D 
and  the  perturbation  h  and  6  /  6 1  is  the  two-dimensional  substantial 
derivation.  V'hen  linearized,  Equation  (6)  becomes 

dt/dt  -  (f  +  Z)  •  db/dt  +  vL  3  H-1(f  +  Z)/3y  =  0  (7) 

where  d/dt  =  31+  U/9x,  Then  introduci::^  the  gecstrophic  approximations: 

u  =  -g'f“18h/3y,  v  =  g'f”^3h/9y.  r,  =  g'f*  V^h  (8) 

and  denoting 

h  =  4>  (y\  expi  a  (jc-ct),  c  =  a  +  ib  (9) 

p  (y)  -  (f  +  7)  Id  (io) 


v/e  then  obtain  from  Equation  (7)  the  eigen  value  equation 


9  "(y)  -  ia  2  +  f  P(y)  /  g'j  9  -  D  P'  9  /  {  TJ  -  c  )  =  0 


(11) 


where  a  prime,  except  to  g,  denotes  differentiation  with  respect  to  y.  The 
boundary  condition  is  that  9  vanishes  at  y  =  ±  »'  •  .  Multiplying  (11)  by  the 
conjugate  eigen  function  4’  *  ancl  averaging  the  result  over  all  values  of  y, 
we  obtain 


J 


DP'(y)  9  9  * 


(U  -  a)2  +  b2 


-1 


dy  =  0 


(12) 


by  equating  the  imaginary  part  to  zero  (Lin,  1955).  Equation  (11)  indicates 
that  if  there  are  unstable  modes  (i.e.  b  >  0),  then  P'(y)  must  change  sign 
within  the  jet.  Therefore,  the  potential  vorticity  of  the  mean  flow  P  (y) 
has  an  extrema  within  the  jet  as  a  necessary  condition  for  instability  of  the  jet. 
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'  here  C,  <<  f  or  the  Rossby  number  is  extremely  small,  we  have 

D  P’(y)  -  U  (U"  /  U  -  fZ  /  g*  D) 


(12 


If  U"/U  is  always  finite,  then  P'(y)  never  vanishes  or  the  jet  is  stable  whe 


fZ.'r'Di  >  cr  c.x.  (U"/  U) 


(14 


where  Dj  -  D  ( o-.  f  ar.d  max. 
jet  v/> ' h  a.  given  velocity  profile 
ottis-  things  bfcir.v  ecu  a;.. 


means  the  maximum  value.  Therefore,  a 
is  the  more  stable,  the  smaller  the  depth, 


if 


V  hen  D  is  sufficiently  large  (F<<  1),  P(y)  0  and  D  P' --.'dZ/ 
=  U"(y).  It  is  easily  see:,  that  Equation  (11)  reduces  to  the  classical  eigen 
value  equation  of  c  two  dimensional  inviscid  disturbance  for  the  mean  flow 
U(y)  v/hich  is  unstable  in  the  vicinity  of  certain  discrete  wave  numbers  (a) 
from  T cllmien' s  theory  l  in,  1 9 5 5) 0  Therefore,  there  must  exist  a  maximum 
crfFc:’  depth  D  j  r  above  which  the  jet  is  unstable  and  the  condition  for  this 
number  is 

0  <  f2  /  (g1  Dlc)  <  max  U"/  U  (1 

The  more  precise  criterion  requires  the  solution  of  Equation  (11) 
which  j.s  reduced  to  the  following  dimensionless  form  (if  R  <<  1): 

$  (S,  )  -  (aD2  +  F)  *  -  (  W"(  4  )  -  W  (t)  j  (V-c0)'S  =  0  (11 

where  I  J 

4  =  y/  L,  U(y)=UQW(e).  cq  =  c/ Uq,  a  q=  a  L  (16 

In  this  equa.tion,  F  is  constant  to  the  order  of  R,  since  the  expansion  of  the 

undisturbed  depth  D  (y»  is  given  by  L  (y)  =  Ec  +  O  (£,  R)  from  (IV-38).  Whe: 

the  mean  flow  is  given  by  IV  =  sech  |  (Bickley  jet),  the  solution  of  (11a) 

2  2 

becomes  <{>  =  sech  £  provided  aQ  =  6  cq  =  (4-F),  Then  long  quasi- 

ionary  waves  (aQ  ~  0,  cu  ~  0)  become  neutral  at  F  ~  4  and  thus, 
unstable  solutions  exist  near  the  neutral  waves  from  Tollmien's  theory 


!• 
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(Lin,  1955),  A  sufficient  condition  for  the  instability  of  this  jet  is  that 

P  <  4.  When  the  physical  constants  are  taken  as  f  =  10  sec  ,  g'  = 

Z  —  1 

1  cm  '  sec  and  -  SO  km,  Dq  =  620  m  corresponds  to  F  =  4.  When 
the  mean  depth  y.x  :ec  ds  t.ais  value,  the  jet  will  be  unstable  provided  R<  <  1 
or  U  <<  5  \ -  "  L.)  m/sec,  S:ern  suggested  that  the  initiation  of  the  Gulf 
Str  .'•••* :  uea.vit  rs  r.  >/  ih  of  Cape  Hatteras  may  be  due  to  the  increase  of 

its  me  a  a  dp.  hi. 

Pedlosky  (1963)  discussed  the  stability  of  a  two-layer  system 
in  which  the  rr.om.i  th.w  ?.r.  the  lower  layer  U~,(y)  is  finite.  Quantities  in 

w 

the  upper  and  iow  :  v  j  <  are  distinguished  with  suffices  1  and  2,  respec¬ 
tively.  The  con-.e-.rv  i; Ion  of  potential  vorticity  and  the  conservation  of  vol¬ 
ume  ar*3  Imcp  c ired  Li  eih  layer  by  the  perturbation  method  as  before,  under 
the  c  o  •  L  ton  a  cf  .re  hydrostatic  pressure  and  geostrophic  balance.  The 
ncr. Vrr  ensicna]  perturbation  equations  for  the  pressure  p^  and  p^  can 
be  written  as 

(d2/d|  2  -  a2)Pi  +  S;  Pi  (Vi  -  c)-l+Fj,  (prp2)(-l)i=0  (i=l,2)  (1 

where  all  the  quantities  are  made  dimensionless  by  scaling  the  variables  with 
(16)  (and  dropping  the  suffix  o  in  a  and  c),  F-  is  the  Froude  number 
defined  by  (lb)  with  Dq  replaced  by  or  the  undisturbed  thickness  of 
the  upper  and  lower  layer,  and 

Si  =  PR"1  -  V^"  (|  )  -  F.  H.'  (4),  (i  =  1,  2)  (18) 

in  which  p  is  the  non-dimensional  rate  of  change  of  f  with  latitudes  and  a 
prime  indicates  the  differentiation  with  £,  ,  Multiplying  each  equation  of 
‘  ■  7)  by  (complex  conjugate  of  p^)  and  subtracting  from  this  its  conjugated 
equation,  we  have 


. 


-32- 


(Pi*Pi*  ~  P^P^))'  +  2  b  S-jp.j  |W.  -  c| 
=  F.(-l)1  (p1P2*  -  P2Pl*)»  (  i  =  1,  2 ) 


-2 


(19 


in  whi~h  b  is  the  imaginary  part  of  non-dimensional  c(=  a  +  ib).  Integrating 
two  equations  of  (18)  between  £  =  -1  and  1,  where  pj  and  p^  vanish  (the 

values  ±1  may  be  replaced  by  any  other  constants)  and  summing  the  results, 


we  have 


„  1  2 


2b  , 

J 


D. 


— 1  i=l 


iSi  |Pi|  Vi‘cl  <U  =0 


(20) 


Therefore,  if  b  ^  0,  S.  must  be  positive  in  some  regions  and  negative  in 
others. 

One  interesting  application  of  Equation  (20)  is  the  case  where 
=  0  and  Vr ^  =  0,  neglecting  p.  Then  in  the  upper  layer 


W"  (£)  -  Fx  Wj  (|)  =  0  (21) 

since  D j '  ( £ )  =  -  W  (£  )  in  non-dimensional  units.  If  there  is  a  wall  at  £  =  0, 
the  mean  flow  is  given  by 

W j  =  exp  £  n/  F j  (~CH<  £  <  °)  (22) 

Since  D2'  (£  )  =  -  Dj'  (£  )  =  Wj  and  S£  does  not  change  sig  n,  Equation 
(19)  yields  b  =  0  and  the  flow  is  stable.  Equation  (22)  represents  the 
velocity  profile  of  the  Gulf  Stream  obtained  by  Stommel  (1958)  with  a 
condition  of  conservation  of  potential  vorticity.  However,  if  we  change  the 
at  bottom  to  a  variable  slope  of  suitable  forms,  S^  may  be  made  both 
positive  and  negative  and  b  may  become  finite.  Therefore,  the  effect  of 
a  variable  bottom  slope  may  change  the  stability  characteristics  of  the  flow 
without  changing  the  kinematics  of  the  mean  flow. 


•  H  !  .  *  '  -  (  , 
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Howard  (1961)  proposed  a  method  of  obtaining  criterion  for 


stability  of  stratified  two-dimensional  shear  flows.  This  method  was  applied 
to  the  jet  in  a  two  -  layered  system  by  Pedlosky  (1963,  1964a),  In  the  follow¬ 
ing  process,  p  will  be  neglected  for  simplicity.  Introducing  the  functions 

=  Pi  (  W.  -  c)"1  (  i  =  I,  2)  (23) 

into  Equation  (19),  multiplying  each  equation  by  integrating  between  -1 


and  1  with  £  and  adding  the  resultant  equations,  we  have 

7^7  -  c)2o7  +  <  v;7-  c)z~z  =  rwr^'/j 

where  a  bar  indicates  integration  between  -1  and  1  and 

2 


(24) 


Qj  F.  =  |  ^i'(  1)1  2  + 

^i  -  j 2  / 2 


j  = 


(25) 

(26) 


Therefore,  and  J  are  positive  finite  quantities.  The  imaginary  and  real 


parts  of  (26)  are  respectively 


(v  !  -  a)  C>!  +  (W2  -a  )  C2  =  0 


(27) 


TV12  +  ^2Z  C2  =  (a2  +  b2)(Q1  +  Q2)  +  (  W2  -  J  (28) 

When  n^  <:  V  •  <:  m^,  the  following  inequality  holds 


- 2  ’ 

T*’  \  °  T 


£  (Wi  -  raj)  -  n.)  Qj  .<  0  (29) 

Integrating  Equation  (29)  between  4  s-1  and  1*  and  substituting  relations 


(26)  and  (27),  we  have 


,  (a^  +  b^)  -  a  (m^  +  n^)  +  mj  n^ 


(QX  +  02) 


+  (mj  +  nj  -  m2  -  n^)  V!  Q2  +  (m2  n2  -  mj  n^)  C2  <  0 


(30) 


in  which  (V;2  -  Wj)2  J  is  dropped  from  the  righthand  side,  since  it  is  positive, 


. 

>• 


in  (30)  ma  be  replaced 
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'"hen  mj  +  +  n2  and  m^  +  nj  >0,  v  2 

by  n2  and  the  inequality  (30)  becomes 

i  a  -  (m,  +  n,)/2  j  +  b^  (mi -nip/ 4 

L  1  J 


+  (mj -n2  )  (nj -n2  )  C2  /  (  Qj  +  Q2  ) 


(31) 


This  equation  indicates  that  the  complex  phase  speed  lies  within  a  se  janL-circle 
as  indicated  in  Fig.  7,  in  which  it  is  assumed  that  m^  »  m2  and  n^  <<  n2. 
It  implies  that  the  flow  is  capable  of  greater  growth  rates  if  the  radius  of 
the  semi-circle  is  larger  (for  instance,  coinparison  of  the  case  n^  >  n2  with 
the  case  n^  <  n^). 

Pedlosky  (1963,  1964b)  also  showed  that  the  theorem  of  Tollmien 
(Lin,  1955)  is  valid  for  a  two-layer  model.  It  is  assumed  that  the  velocity  of 
the  lower  layer  is  zero  and  cg  is  the  value  of  VI  j  at  the  point  where  Sj  =  0 
or  the  vorticity  has  an  extremum  with  =  0.  The  K-  functions  of  Lin's 
method  (Lin,  1955)  for  the  two  layers  are  defined  by 

Kj  (|)  =  Sj  (V  .  -  c g) 1  ,  (i  =  1,2)  (32) 


where  W 2  =  0.  If 

Kj  +  K2  >  tt2  /  (  2  a2)  (33) 

where  2a  is  the  width  of  the  jet,  it  is  shown  that  Equation  (16)  has  a  neutral 
solution  with  phase  speed  cg.  The  proof  is  given  by  the  variational  method 
fcr  the  functional  M  defined  by 


:r 


-1 


-1 


M  =  1  *3  KilpjTFj  -  Opj/as^Fj  -  tm-Ffc)  ^Pi  Ft  •  (3- 

<_  J 

where  a  bar  denotes  the  integral  with  ^  between  -1  and  1  as  before  and 


the  summation  is  for  i  =  1  and  i  =  2.  The  functional  Vl  has  a  maximum 


.  ■  •  t  . 
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where  its  variations  with  and  6  .A  I  5  p-  and  6  M/  6  p^  vanish.  If 

c  =  c  ,  these  variational  conditions  are  equivalent  to  Equation  (17).  There- 
s 

2 

fore,  a  neutral  solution  p^  g  exists  with  an  eigen  value  a  g  =  max.  vl  for 
c  =  cct  It  is  also  shown  that  the  flow  becomes  unstable  to  wave  numbers 
slightly  less  than  a  by  use  of  a  method  modified  from  the  one  developed  by 
Lin  (1955).  Let  p^  be  an  eigen  function  for  a  which  is  near  a  s.  Multiply¬ 
ing  Equation  (17)  for  p^  g  and  p^,  respectively  with  p^  and  p^  g  )  summing  all 
of  them  and  integrating  with  6.  over  the  range  (  -1,  1  )  we  have 


a  ^  a 


')/( 


c  - 


c  )  = 
s  ' 


- - n - Ti  /  - 

sVipisFi  <  Wi"c  >  /  2  Pipi 


sFi 


-T 


(35) 


By  taking  limit 


2  2 

a  — >  a  „  ,  c  — ;c0,  we  can  show  that 

*  S  o 


(da  */  dc)  =  A  +  iB 

c  — i  cs 


(36) 


with  E  >  0.  Thus,  if  a  2  decreases  slightly,  the  imaginary  part  of  c  becomes 


“1 


po 


sitive.  It  is  shown  that  the  sine  type  jet  !  {Vf  cos  (  it  ^  /  2  a)  '  or  the 


Gaussian  jet  !  W  '  exp  (  -  r^  ) 


has  a  neutral  solution  and  the  unstable 


regime  near  it. 

Theories  of  stability  of  an  atmospheric  zonal  flow  were  reviewed 
by  Phillips  (1963)  for  geostrophic  disturbances  and  also  by  Miles  (1964)  for 
disturbances  of  more  general  types.  The  latter  showed  that  disturbances  of 
sufficiently  short  wavelength  in  flows  with  typical  velocity  distributions  are 
always  unstable  but  that  they  have  a  rate  of  growth  vanishing  directly  with 
the  wave  length.  Also,  the  latter  proved  that  flows  with  sufficiently  small 
velocities  are  always  unstable  for  all  wave  lengths. 
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6,  Rossby  v  aves  as  a  Cause  of  Meanders 

Charts  of  surface  isotherms  of  the  western  North  Pacific  for 
each  ten-day  period  have  been  made  by  the  Japan  Meteorological  Agency 
since  about  1950  by  use  of  the  data  collected  by  merchant  ships  and  fishing 
vessels.  Comparison  of  the  surface  isotherms  with  the  current  pattern 
based  on  dynamic  topography  and  direct  measurements  with  G  E  K  showed 
that  the  strong  gradient  of  the  isotherms  almost  corresponds  to  the  path  of 
the  Kuroshio  particularly  in  the  area  east  of  Japan  where  the  Kuroshio  flows 
eastward  after  leaving  the  coast.  In  many  instances,  these  isotherms  show 
meanders  of  wave  lengths  of  a  few  hundred  to  more  than  a  thousand  nautical 
miles  with  amplitudes  up  to  300  miles  (Ichiye,  1957).  In  most  cases,  the 
wave  pattern  of  meanders  seems  to  move  eastward  in  the  direction  of  the 
mean  flow  although  the  speed  seems  to  be  much  less  than  the  mean  current. 
Also,  the  meanders  of  the  isotherms  usually  showed  a  feature  of  tilting 
troughs  and  ridges  in  the  direction  from  southwest  to  northeast.  The 
harmonic  analysis  in  space  of  the  15°C  isotherm  indicates  that  a  large  ampli¬ 
tude  of  meanders  with  a  wave  length  of  1000  n,  miles  breaks  into  waves  of 
higher  harmonics  after  about  one  month  or  three  consecutive  charts.  These 
features  suggest  that  the  meanders  of  the  surface  isotherms  are  quite  similar 
to  meanders  of  the  westerlies  in  the  upper  atmosphere  indicated  by  daily 
pressure  charts,  Rossby  (1939)  introduced  an  idea  of  planetary  waves  to 
explain  the  dynamics  of  westerlies  meandering  and  many  meteorological 
works  were  devoted  to  development  of  his  theory  since  then,  Ichiye  (1950, 
1956,  1957)  tried  to  explain  various  features  of  the  meanders  of  isotherms 
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from  the  planetary  (Rossby)  wave  theories,  for  instance,  the  breaking  of 
meanders  into  higher  harmonics  from  the  dispersion  theory  of  Yeh  (1949)  and 
the  tilting  of  troughs  from  Viachta's  model  (1949),  However,  the  original 
model  of  Rossby  waves  seems  to  be  inadequate  for  the  meanders  of  the 
western  boundary  currents  even  when  they  become  zonal  flows  after  leaving 
the  continent,  because  the  effect  of  variation  of  relative  vorticity  within  the 
current  or  branch  of  multiple  currents  is  much  larger  than  p-effect,  There¬ 
fore,  in  most  cases  meanders  of  the  Gulf  Stream  may  be  treated  as  perturba¬ 
tions  of  a  shear  flow  in  a  system  uniformly  rotating  around  a  vertical  axis, 
as  discussed  in  the  preceding  sections. 

On  the  other  hand,  when  we  consider  the  Gulf  Stream  Extension 
as  a  portion  of  the  general  circulation  of  the  North  Atlantic,  the  p-effect 
cannot  be  neglected.  Fofonoff  (1961)  showed  that  an  eastward  mean  flow  is 
accompanied  with  wavy  perturbations  due  to  stationary  Rossby  waves,  but  for 
a  westward  flow  the  perturbations  are  not  of  an  oscillatory  type.  He  suggested 
that  in  a  model  of  oceanic  circulation  of  a  rectangular  basin,  a  basic  flow  in 
its  southern  half  is  almost  zonal  and  westward,  while  a  basic  eastward  flow 
in  its  northern  half  is  superposed  with  stationary  Rossby  waves.  .Toore  (1963) 
discussed  the  stationary  Rossby  waves  in  the  interior  part  of  the  ocean 
including  the  effect  of  eddy  viscosity.  Both  authors  treated  the  barotropic 
model. 

First,  Fofonoff's  (1961)  theory  is  discussed.  In  a  homogeneous 
ocean  of  a  constant  depth,  components  of  two-dimensional  stationary  motion 
can  be  defined  by  a  stream  function  such  that 
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u  =  9  ^  /  9  y,  u  =  -  9  ijj/  3x  (1  a,  b) 

Throughout  this  section,  the  x-axis  is  taken  eastward  and  the  y-axis 
northward.  The  equation  of  absolute  vorticity 

Z  =  I,'  +  f  =  f  -  V  (2) 

is  given  by 

u  9  Z  /  8  x  +  v3Z./9y  =  0  (3) 

-where  t,  is  the  relative  vorticity.  Substitution  of  (1)  into  (3)  leads  to 

94j/9y*  aZ/9x-94>/9x‘  9Z/9y  =  0  (4) 

Therefore,  Z  becomes  a  function  of  In  the  simple  case  of  intertial 
flows,  Z  is  assumed  to  be  a  linear  function  of  Further  f  is  assumed 
to  be  fQ  +  p  y.  Then,  we  have 

f0  +  py  -  V2+  =  co  +  c,4<  (5) 

where  CQ  and  are  constants.  The  constant  Cj  is  determined  from  a 
condition  that  the  velocity  is  zonal  and  equal  to 

9  v{j  /  9  y  =  U  =  p/Cj  (6) 

wherever  l  vanishes.  The  constant  Cc  is  taken  equal  to  fQ,  Then 
Equation  (5)  becomes 

(U/  (3)  V24)  +  4<  =  Uy  (7) 

The  ocean  is  assumed  to  be  bounded  by  x  =  ±  Lj  and  y  =  0,  L.  The  boundary 
condition  is  that  i|i  =  0  along  the  boundary.  Introducing  non-dimensional 
quantities 

x'  =  tt  x/  L,  y'  =  uy/  L, 

4>'  =  /  |  ul  U  kG  =  (1/tt)  (h  L2/  iul)1/2  »  1  (8) 

and  dropping  primes,  we  obtain  the  non-dimensional  form  of  (7) 
i  **  3  2 

kQ  V  \\>  ±  4>  =  y 


(9) 


. 

. 
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where  the  plus  sign  is  associated  with  an  eastward  flow  (positive  U)  and  the 
minus  sign  with  a  westward  flow.  The  boundary  condition  is  that  ^  =  0  at 
y  =  0,  tt  and  at  x  =  ±  a  =  *  tt  Lj  /  L. 

It  is  instructive  to  discuss  the  solution  of  (9)  in  relation  to  the 
Rossby  waves.  The  homogeneous  solutions  of  Equation  (9)  are 

Q  =  Aq  sinh kQ  y,  =  An  cosh  kn  x  sin  ny  (10  a,  b 

for  the  westward  flow,  where  n  is  integers  and 

2  2  7 

kn  =  k^  +  n  (11) 

with  n  'Z,  1,  The  solution  of  (9)  satisfying  the  boundary  condition  can  be 


obtained  as  follows:  Aq  is  determined  so  that  vanish  at  y  =  0  and 

tt,  where  i^p  =  **  Y  is  particular  solution  .  A^  is  determined  by  a 
Fourier  expansion  of  <\iQ  +  at  x  =  a.  The  complete  solution  is  given  by 
tp  =  -  y  +  ir  sinh  kQ  y/sinh  kQ  it 


+  2  C  ("l)n+^  kA^  cosh  k  x  sinny 
'  ’  o  n  3 

n=  1 


n  (kQ^-  n2)  coshk^a  ' 


w  <X  { 

j 


(12) 


The  solution  for  an  eastward  flow  (positive  sign)  is 
4/  =  y  -  tt  sin  kQy  /  sin  kQ  it 


+  2  E  (-l)n+*k  ^  cos  k  x  sin  ny  n  (k  ^  -  n^)  cos  k  a 
'  ’  o  n  J  o  n 

n=  1  ’ 


-1 


where 


k  2  =  k  2  -  n2 
n  o 


(13) 


(14) 


Equation  (12)  becomes  simpler,  because  k  is  large  compared  with  unity 
for  the  reasonable  set  of  constants.  For  instance,  kQ  -  15,9  for  L  =  2500  km, 
=  2  x  10"13  (cm  sec)“l  and  U  =  -  5  cm  /  sec.  Therefore,  coshknx  =  cosh  ko  x 
for  n  <  k0,  and  the  summation  term  of  Equation  (12)  is  replaced  by 


• '  J 

, 


.  ,  J  =  V 
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(  4jq  +  )  x  cosh  kQx/cosh  kQa.  Further  approximation  for  hyperbolic 


functions  leads  to 


tjj  ~  -  y  -  tt  exp 


<  ) 

-  exp  -  kQ  (a  +  x  )  >  j 


(15) 


This  solution  gives  the  circulation  pattern,  which  has  a  narrow  eastward  jet 
with  a  width  of  L/  it  kQ  (  -  50  km)  along  the  northern  boundary  and  a  broad 
westward  current  in  the  whole  interior  area.  (Fofonoff,  1954).  Solution  (13) 
for  positive  U  suggests  that  the  streamlines  consist  of  a  complex  series 
of  eddies  and  that  no  eastward  jet  is  evident.  This  solution  can  be  inter¬ 
preted  as  superposition  of  stationary  Rossby  waves  on  the  mean  westward  flow. 


A  brief  discussion  of  the  barotropic  Rossby  waves  will  be  presented 


here.  The  linearized  two-dimensional  momentum  and  continuity  equations 
of  homogeneous  water  with  a  uniform  zonal  current  U  (x-  direction)  are 
respectively  given  by 


(16) 


(  a/3  t  +  U  8  /  9x)  n  +  H  div  “v  =  0 


(17) 


/  / 

where  v  and  k  are  velocity  vector  and  unit  vector  of  the  vertical  direction; 


t-  is  the  surface  elevation  and  H  is  the  mean  depth  assumed  as  constant. 


A  dimensionless  notation  (primed)  is  introduced  in  Equations  (16)  and  (17) 


in  terms  of  characteristic  parameters: 


(x,  y)  =  (x',  y')  L,  rt  =  'H 


(u,  v,  )  =  (u1,  v')  U,  t  =  (L/U)f 


P  =  8f  /  9y  =  p*  (L/  f) 


(18) 


' 

t 
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When  a  prime  is  dropped,  Equations  (16)  and  (17)  can  be  written  as 

R  6  <"  /  6  t  -  (l  +  py)Tc'xv>=  ~  G  V  q 

6  q  /  6  t  +  div  v  =  0 


where 


R  =  U  /  f  L,  G  =  gH/(fUL) 
6/  6t  =  9/  at  ±  9  /  8  x 


(19) 

<20) 


(21a,  b) 
(22) 


in  which  a  positive  sign  corresponds  to  an  eastward  flow  and  a  negative  sign 
a  westv/ard  flow.  The  numerical  constants  used  in  the  inertial  circulation 
(  L  =  2500  km,  f  =  10** 4  sec  U  =  5  cm/sec,  p  =  2  x  10  ^3  (cm)  *)  with 
H  =  4000  m  give 

R  =  2  x  10”4,  G  =  3  x  103,  p  =  0.  5  (23) 

- 'y 

The  quasi-geostrophic  relation  for  v  to  the  order  of  R  is  given  by 

(1  +  Py)~v  /  G  =  k  x  V  r|  -  R  (1  +  py)**1  (  6/  6t)Vq  (24) 

Substituting  (2  3)  into  (19)  we  have 

R  (  6/  6t)  (  V2q  "2  hy  /RG)  +  Pqx  =  0  (25) 

in  which  the  term  1  +  py  is  replaced  by  1.  When  q  is  assumed  to  be  of 
a  form  expi  (  <rt  +  kx  +  £  y).  Equation  (24)  yields  the  frequency  equation 

R  w  (k2  +  £  2  +  2  pH  i  +  M)  -  pk  =  0  (2  6) 

where 

M  =  (  RG  )~  1  to  ss  cr  ±  k  (27  a,  b) 

In  Equation  (27b)a  plus  sign  corresponds  to  an  eastward  flow  and  a  minus 
sign  a  westward  flow.  The  waves  which  satisfy  the  frequency  Equation  (26) 
are  called  Rossby  waves.  Equation  (26)  indicates  that  if  the  waves  are  stable, 
the  amplitude  should  be  changed  exponentially  unless  the  mode  is  purely  zonal 
(  f  =  0)  .  V  hen  k  is  assumed  real  and  £  =  v  -  i  k  »  Equation  (26)  yields 
R  to  ( k3  +  +  P^  +  ivl)  -  pk  =  0.  a  -  p 


(28  a,  b) 
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The  waves  become  stationary  when  <r  =  0.  However,  the  waves  with  a 
constant  amplitude  are  possible  only  for  the  eastward  flow  (positive  sign 
in  (27  b).  In  this  case  (28a  )  becomes 

k2  +  vl  -  p/  R  (29) 

by  neglecting  terms  of  the  smaller  order.  Equation  (29)  is  equivalent  to  the 
characteristic  Equation  (14),  since  pR"  is  equal  to  p  L.‘J  /  U  (  =  ir  ),  On 
the  other  hand,  when  cr  =  0  for  the  westward  flow,  Equation  (28a)  is  satis¬ 
fied  only  when  k  is  imaginary.  If  we  put  k  =  i  k1,  (27a)  becomes 

VZ  +  p  /  R  -  (k')2  (30) 

which  is  equivalent  to  (11),  It  is  easily  seen  that  in  this  case  the  eigen  func¬ 
tions  have  an  exponential  function  of  x  as  a  factor, 

Toore  (1963)  treated  a  model  of  stationary  wind-driven  ocean 
circulation  in  a  two-dimentional  homogeneous  system  including  the  inertia  and 
eddy  viscosity  terms.  The  velocity  can  be  expressed  by  relations  (la,  b) 
using  a  stream  function  ip  under  an  assumption  of  non-divergent  flow’.  Since 


8(fip)/9y  =  fu+p4'  (31),) 

the  Navier-Stokes  equations  in  a  p-plane  can  be  written  as 

u3u  /3  x  +  v9u/  9y=  -  p  "^p'/Sx  +  W  +  K  V2  u  (32) 

u3v/3x+v9v/3y  -  pip  =  ~p"'^3pl/3y  +  K'v'  v  (33) 

by  introducing  a  modified  pressure  p1  defined  by 

p'  =  p  +  p  f  vp  (34) 


In  Equation  (32)  W  is  a  zonal  wind  stressfunction  expressed  as  a  body  force. 
Cross  differentiation  of  (32)  and  (33)  yields  the  vorticity  equation 

c  (9  4j  /  3  y)  9  ( V2  )/  9  x  -  e  (  3  4;  / 3  xJ»V2  *P)  /  3  y  +  3  ip  /  3  x 
=  d  Vv  /  dy+6V4ip 


(35) 
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after  scaling  the  variables  by  use  of  (8)  and  dropping  a  prime.  In  this 
equation,  new  parameters 

€~U/(j3jL2)  6  =K/(PL3)  (36  a,  b) 

are  introduced  (  e  /  6  is  the  Reynolds'  number  in  a  conventional  sense  ) 
and  the  non-dimensional  wind  function  W  (with  a  prime)  is  equal  to 
V7  /  (LpU),  The  wind  stress  torque  dT?/V  dy  has  a  large  negative  value  in 
an  area  between  20°  N  and  35°  N  in  which  the  ocean  circulation  has  a  rather 
meridional  structures  In  the  north  and  south  of  this  area,  the  current  shows 
a  strong  zonal  structure  and  the  wind  stress  torque  is  small.  Therefore,  we 
can  neglect  this  term  for  the  flow  to  be  discussed.  This  assumption  is  contra¬ 
dictory  to  a  model  of  wind-driven  circulation  treated  by  Carrier  and  Robinson 
(1962)  whose  solution  shows  a  strong  and  narrow  inertial  jet  at  the  latitude  of 
maximum  wind  stress  torque  dividing  anti-cyclonic  gyre  to  the  south  and 
cyclonic  gyre  to  the  north.  However,  the  problem  to  be  discussed  here  is 
stability  of  a  zonal  current,.  Further,  it  is  assumed  that  the  zonal  flow  has  a 
feature  of  a  boundary  layer  current,  that  is, 

=  0,  d2^;  /  9  y2  =  3  u  /  8  y  =  0  on  y  =  0  (37  a,  b' 

and  on  y  =  y^  (  <  1  )  becomes  the  solution  of  (35)  in  the  interior  ocean 
in  which  the  inertia  and  viscosity  terms  may  be  neglected,  but  the  wind  stress 
torque  is  finite.  Boundary  conditions  (37  a,  b)  indicate  that  the  flow  is  purely 
zonal  on  y  =  0  and  neither  mass  nor  momentum  is  transferred  across  y  =  0. 
(This  boundary  is  not  necessarily  a  solid  barrier  but  may  be  taken  as  the 
center  of  the  jet).  From  conditions  (37  a,  b),  the  stream  function  is  assumed 
to  be  of  a  form 

4»  =  y  F(x)  +  y3  F3  (x)  +  ...  (38) 
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Substituting  (38)  into  (35)  and  equating  the  coefficients  of  y  to  zero,  we  have 
c  (F  F"'  -  F*  F*»)  +  F*  =  6FM"  (39) 

where  a  prime  denotes  differentiation  with  x.  Once  integrating  (39)  with  x, 
we  have  p 

6  i  F  F"  -  (  F'  )2;  +  F  =  -  C  +  6  F"«  (40) 

where  C  is  an  integration  constant.  The  solution  of  (40)  is  given  by 

F  =  -  C  +  C  exp  (  -  mx)  (41) 

provided  that 

m3  -  6  e  C  m2  +  6=0  (42) 

V7hen  6  <<  1,  Equation  (42)  is  reduced  to 

€  Cm2  -  1  =  0  (43) 


Equation  (43)  is  the  same  as  Fofonoff’ s  (1954)  characteristic  equation  for  an 
inertial  circulation,  equivalent  to  (11)  and  (14).  V7hen  C  >  0,  m  is  real  and 
Equation  (41)  indicates  a  westward  flow  which  vanishes  at  x  =  0,  decreasing 
uniformly  from  further  east  (  x  =00)  if  m  is  taken  positive.  "V:hen  C  <  0, 
a  mean  flow  is  eastward  but  m  is  imaginary  and  the  flow  is  oscillatory. 

General  Equation  (42)  has  real  roots  only  if 

4  6  “2  (  e  C)3>  27  (44) 

whilst,  in  the  contrary  case,  it  has  a  pair  of  conjugate  complex  roots.  It  is 
easily  seen  that  the  westward  flow  (  C>  0)  always  has  complex  roots.  Since 
we  may  take  |  C  |  as  unity,  relation  (44)  becomes 

4  U3  >27  PK2  (45) 

for  the  westward  flow.  It  is  noted  that  the  scale  distance  L  is  not  included 
in  this  criterion.  For  U  =  10  cm  /  sec,  this  inequality  is  valid  if  K  <  2.7  x  10^ 
cm^/ sec,  A  simple  algebraic  calculation  shows  that  the  roots  of  (42)  consist 
of  one  negative  and  two  positive  values  if  (45)  is  valid  and  one  negative  and 
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complex  conjugate  roots  with  a  positive  real  part.  However,  the  negative 
root  for  an  eastward  current  is  physically  possible  only  when  the  number  6 
becomes  of  the  order  of  unity  or  more  or  the  eddy  viscosity  is  sufficiently 
large,  because  otherwise,  the  absolute  value  of  the  negative  root  is  so  large 
that  the  velocity  becomes  unreasonably  large,  as  x  decreases  from  zero. 

For  instance,  for  L  =  100  km  and  U  =  10  cm/ sec,  6  =0,05  and  0,5 
respectively  for  K  =  10^  and  10^  cm^/sec,  and  the  negative  root  for  an 
eastward  flow  is  -6.  1  and  -2,4,  respectively.  On  the  other  hand,  the  nega¬ 
tive  root  for  a  westward  flow  is  -2,0  and  -1,7  for  the  respective  value  of 
K,  Therefore,  for  the  westward  flow  the  stationary  current  decreasing  from 
east  to  west  is  possible  for  a  relatively  low  eddy  viscosity, 

Moore  (1963)  further  solved  Equation  (35)  for  the  wind  stress 
function  W  •-  -  cos  tt  y  /  L' .  He  linearized  the  equation  in  a  way  similar  to 

the  boundary  layer  theory,  assuming  that  the  main  interior  current  flows  in 
the  zonal  direction,  thus  obtaining 

-  €  U  ( y)  ^ 1 1 1  +  ^ '  =  d  V.7  /  dy  +  6  ^  M 1 1  (46) 


where  U  (y)  is  a  non-dimensional  zonal  velocity  and  is  parameterized 
though  it  depends  on  y.  For  the  realistic  constants  UQ  =  10  cm  /sec,  p  =  10“*- 

(cm  sec)“*,  K  =  2  x  10^  cm2/ sec  and  L  =  5000  km,  corresponding  to 

- 1  2 

6  c  =  2,5  x  10  ,  he  obtained  the  streamlines  showing  damped  Rossby 
waves  in  the  western  part  of  the  northern  half  of  the  ocean,  suggesting  that 
the  observed  meanders  of  the  Gulf  Stream  rather  may  be  a  feature  of  the  stable 


solution  than  due  to  instability. 


Il'in  and  Kamenkovich  (1963)  discussed  Equation  (35)  by  setting 
3/2 


6  =  o  € 


(47) 
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where  e  is  a  positive  constant,  since  in  general,  e  >>  6*  t  urther,  ii  the 


wind  stress  W  is  a  linear  function  of  y,  it  is  assumed  that 

dW/dy  =  b;  4>=W(y)F(x)  (48) 

Then,  Equation  (35)  yields 

-e  e  3^2  F""  +  6  b  (  F'  F"-  FF1"  )  +  F'  +  1  =0  (49) 

instead  of  Equation  (40).  When  an  ocean  is  bounded  by  x  =  0  and  x  =  1, 
the  boundary  conditions  at  these  points  are 

F  =  F'  =  0  (50) 


By  use  of  boundary  layer  techniques,  F  (x)  is  expanded  into  an  asymptotic 
series  by  powers  of  e  like 

F  (x)  =  Fq  (x)  +  F1  (x)  +  e  F^  (x)  +  .  . . 


where 


+  (q)  +  n/c  E^(q)  +  ...  +GQ(r)+\/€Gj(r)  +  ,..  (51) 

q  =  x/n/T,  r  =  (  1  -  x)  l  \J1  (52  a,  b) 


and  En  or  Gn  is  significant  only  near  x  =  0  or  x  =  1,  respectively. 
Substituting  (51)  into  (49)  and  equating  coefficients  of  the  same  power  of  €, 


we  have 


F0*  +  1  =  0  ... 

-eE  +  b  I" E  '  E  ' 1  -  (  F  (o)  +  E  )  E 
o  j  o  o  'o'  o  o 


+  Eq  =  0 


(53) 

(54) 


-e  G  -  b  Fg  1  G^’*  -  (Fn(l)  +  G  )  G  ,n]  -  G'  =  0  (55) 

o  o  o  '  o'  f  o  o  !  o  ' 


where  a  prime  indicates  differentiation  with  x,  q  and  r,  respectively.  The 
boundary  conditions  become 


E 


n(°)  +  Fn(o)  =  0,  En'(o)  =  0,  En+i'  (o)  +  Fn'  (o)  =  0  (56) 


G„(o)  +  Fn(l)  =  0,  Gq'  (o)  =  0,  -  Gn+1’(o)  +  Fn'  (1)  =  0  (57) 
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In  expansion  (51),  En(q)  and  G  (r)  must  decrease  exponentially  as  q 
and  r  become  infinite.  Therefore,  in  the  range  of  large  q  and  r,  non¬ 
linear  terms  of  (54)  and  (55)  are  neglected  and  we  have 

-e  E  " '  -  b  F  (o)  E  "  +  E  =  0  q  »  1  (58) 

o  o'  '  o  o  ^ 

-e  G  " 1  +  b  F  (1)  G  n"  -  G  s  0  r»l  (59) 

o  o  °  ° 

From  Equation  (53),  FQ  (x)  =  1  -  x.  Then  G  ~  0,  ^wing  to  the  conditions 
that  G  (o)  —  G0' (o)  —  0. 

inhen  we  assume  that  GQ(q)  is  of  a  fcrm  exp  ^  q,  y  is  t  vo 
rc^ts  ^f  the  cubic  equation 

+  b)(^-l  =  0  (60) 

with  Re  x  <  0.  7  hen  b  >  0  (corresponding  tu  the  southern  half  of  the 

northern  ocean),  we  have 

Xj  =  -l/^b+  0(e),  X2  =  ”b/e  +  0(e)  (61a,  b) 

f^r  small  values  A  e.  These  r^ts  indicate  that  the  boundary  layer  has 
two  kinds  with  thickness  A  the  -rder  cf  \j  e  and  e  n/~c  ,  respectively.  v  hen 
b  <  0  (northern  half),  ve  have 

Xj  =  i  /  N/“l~b|  -  e  /  2  b2  +  O  (e2),  X?  =  "  ‘  /’'Mbl  -  e/2b2  +  O  <e2) 

(62) 

Lr  small  values  ^.f  e.  The  boundary  layer  is  larger  than  in  the  former  case 
and  of  the  -rder  of  \Tl/e,  but  the  functions  show  damped  •_  scillation  with  a 
variable  q.  This  corresponds  to  the  Rossby  waves  in  an  eastward  current, 
as  discussed  by  F^fon^ff  (1961)  aad  vR  ore  (1963)  f.r  a  pure  inertial  flow.  It 
is  easily  seen  fr^m  Equation  (49)  that  this  model  becomes  Carrier  and  Rcbinsoj 
(1962)  model  for  small  values  .  f  e  and  flunk's  (I960)  model  for  large  values 
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of  e.  Therefore,  Equation  (54)  is  essentially  the  same  as  those  treated 
by  jvloore  (1963). 

Welander  (1961)  treated  baroclinic  Rossby  waves  propagating 

in  a  zonal  ocean  current  decreasing  exponentially  with  depth.  He  used 

the  x-y-p  system  with  the  vertical  velocity  <o  =  3p/8t  and  9  (=g  z) 

the  geopotential  of  a  pressure  surface.  It  is  assumed  that  the  perturbation 

is  geostrophic  although  the  basic  current  is  zonal  and  that  a  hydrostatic 

balance  exists  in  the  vertical  direction.  Then,  equations  of  motion  in  the 

horizontal  and  vertical  directions  are,  respectively. 

-fv  =  -  39  /  9  x,  fu  =  -  9<}  /  3  y  (63) 

a  =  -  /  3  x  ( 64) 

where  a  is  the  specific  volume  (a  =  p  **1).  Equations  of  continuity 

and  conservation  of  mass  are  given  by 

9u/9x+0v/3y+3u)/9p  =  O  (65) 

8a  /  3  t  +  U3a  /  0  x  +  f  U  v  +  A  to  =  0  (66) 

P  P 

All  the  dependent  variables  in  (63)  to  (66)  are  for  the  perturbation  except 
U  and  A,  which  are  respectively,  the  zonal  velocity  and  specific  volume 
of  the  basic  system,  satisfying  the  relations 

fU*-*  ,  A  =  -  4  p  (67) 

Derivatives  of  the  basic  quantities  are  denoted  by  a  subscript.  Equation  (66) 
can  be  derived  as  the  perturbation  of  the  total  equation; 

(9/0t  +  u0/0x  +  v9/0y+u9/9p)a  =  0  (68) 

in  which  dependent  variables  are  the  resultant  of  the  basic  quantities  and 
their  perturbations.  Solving  00  from  (66)  by  use  of  '63a)  and  (64),  we  have 
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u>  =  A  -1  (929  /  3p9t  +  U  /  8x0p  -  U  94/ Ox  )  (69) 

P  * 

Substitution  of  velocities  (63  a,  b)  into  (65)  yields 

9  o>  /  9  p  -  pr234  /  9x  =  0  (70) 

retaining  the  p  -  effect,  Since  cu  is  expressed  in  terms  of  function  4 
with  (69),  (70)  represents  the  perturbation  equation  for  9  ,  The  boundary- 
condition  is  given  by 

to  =  0  at  p  =  0;  4  =  p  at  p  =  ''.jc  (71a, 

where  the  condition  at  p  =  0  is  expressed  in  terms  of  4  by  use  of  (69). 
Introducing  a  wave  form 

=  4  *  (y>  p)  exp  ik  (x  -  c  t)  (72) 

Equation  (70)  and  boundary  conditions  (71  a,  b)  transform  into 

(U-c>8(A  *1a<f*/8p)-  !(U  A  "*)  +  Pf2  j  s>  *  =  0  (73) 

p  ;L  p  p  P  J 

(  U-c  )  9  4  */  9  p  =  11^4  *  at  p  =  0;  *  s  0  at  p  =  ^  (74) 

Two  features  can  be  seen  immediately  from  (73)  and  (74  a,  b).  Since  k  does 
not  occur  in  these  equations,  the  wave  speed  is  independent  of  k.  Also,  the 
dependence  of  U  and  A  on  y  can  be  treated  parametrically,  since  no 
y- derivatives  appear. 

The  basic  quantities  are  assumed  to  be 

U  =  UQ  exp  (-p/  pQ),  A  s  Aq  +  Aj  exp  (-p/pQ)  (75  a,  b’ 

where  pQ  is  the  scale  pressure.  Introducing  non-dirnensional  parameters 

p  =  p/pQ»  c*  =  c  /  U0,  7  =  PAlPo/(f2U0)  (7 6  a,  b,  c) 

(73)  and  (74  a,  b)  change  to 

r  1 

,  exp  (-P)  -  c*  ,  d  (exp  P  d  9  *  /  d  P)  /  dP‘  +  t  f  *  =  0 


(77) 


. 


-50- 


(1-c*)  d4>  *  /  dP  +  =  0  at  P  =  0;  c{;*  =  0  at  P  =  CO  (78  a,  b) 

When  c*  =  0  (stationary  waves),  the  general  solution  has  a 

form  exp  k  P  from  which  the  frequency  equation  is  given  by 

K  (k  +  1)  +  T  =  0  (79) 

Wien  7  <  0  corresponding  to  negative  U  (westward  flow),  one  root  of 

o 

(79)  is  positive  and  cannot  satisfy  (78)  and  the  elementary  solution  for  another 
root  cannot  satisfy  (78  a).  Therefore,  no  stationary  waves  occur.  If 
T  >  0,  both  roots  always  have  negative  real  parts  and  the  stationary  waves 
are  possible.  For  t  >  1/4,  the  vertical  structure  is  of  exponential  type;  for 
T  >  1  /4  of  oscillatory  type.  In  the  real  ocean,  ^  is  small  and  the  former 
type  is  usually  possible. 

In  a  general  case  of  c*  ^  0,  the  new  variable 

s  =  exp  (  -P)/  c*  (80) 

transforms  Equation  (77)  and  (78  a,  b)  into 

s  (  s-1  )  d2  <p  *  /  d  s2  +  t  <J>  *  =  0  (81) 

(l-c*)d4*/ds-c*(p  =  0  at  s  =  1/c*;  *  =  0  at  s=0  : 

(82  a,  b) 

The  solution  satisfying  (82  b)  is  expressed  in  terms  of  the  Legendre  function 
which  diverges  at  s  =  1,  except  when  j  =  «  n(  n+ 1)  ( n=l,  2,  3, , )  and  the 
solution  becomes  a  polynomial.  In  general,  therefore,  c*  >  1  in  order 
to  have  a  finite  solution.  Only  solutions  for  7  <-2  i,e.  for  a  westward 
current  are  found.  At  middle  latitudes,  the  physical  constants  lead  to 
T  ~  0,1,  Therefore,  there  are  no  traveling  waves  in  a  westward  current. 

In  the  polynomial  case  (for  an  eastward  flow)  waves  remain  stationary  rela¬ 
tive  to  the  surface  current.  Therefore,  as  in  a  barotropic  model  discussed 
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before,  the  Rossby  waves  of  stable  nature  cannot  propagate  in  the  westward 
flow  in  a  baroclinic  model  with  the  density  varying  exponentially  with  depth. 
It  will  be  worthwhile  to  develop  a  similar  theory  in  a  two-layer  system. 
Also,  it  is  necessary  to  check  the  effect  of  horizontal  shear  of  the  basic 
current  U(y)  by  including  advective  terms  in  Equations  (63  a,  b). 

Ichiye  (1958)  discussed  a  problem  of  long  waves  including 
Rossby  waves  generated  by  the  zonally-moving  wind  stress  of  a  limited 
meridional  extent  in  a  two-layer  ocean.  There  is  no  basic  current  in  his 
model.  His  result  showed  that  the  waves  become  resonant  (unstable)  for 
a  disturbance  propagating  westward  in  a  range  of  Rossby  waves  while  there 
is  no  resonance  for  an  eastward  disturbance.  It  is  yet  to  be  studied  the 
response  of  an  eastward  or  westward  flow  to  a  moving  disturbance  of  a 
pressure  system  or  a  wind  system. 
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7.  Application  of  a  principle  of  Constant  Absolute  Vorticity  Trajectories. 

All  the  theories  discussed  in  the  previous  sections  deal  with 
the  meanders  as  a  perturbation  of  the  basic  current  and  thus  are  based 
on  the  assumption  that  the  amplitude  of  meanders  are  small  enough  to 
neglect  its  square  or  higher  powers.  However,  the  observations  indi¬ 
cated  that  meanders  of  the  Gulf  Stream  and  the  Kuroshio  seemed  to  be, 
on  many  occasions,  much  larger  than  assumed  in  small  amplitude  theories, 
Rossby  (1940)  proposed  a  theory  on  meanders  of  a  narrow  planetary  jet 
from  conservation  of  the  absolute  vorticity.  His  theory  was  applied  to 
forecasting  change  of  the  stratospheric  jet  stream  in  meteorology  and 
was  known  as  the  C  A  V  method,  Ichiye  (1955b)  applied  this  method 
to  courses  of  the  Kuroshio  to  the  east  of  142°  E  in  five  surveys  of 
different  seasons  in  1953  to  1954.  He  used  the  velocity  distributions 
determined  with  dynamic  calculations  based  on  the  800-decibar  level  to 
compute  the  volume  and  momentum  transport  of  the  Kuroshio.  His 
results  showed  that  the  paths  of  the  current  were  within  the  first 
half  wave  length  (300  to  500  km)  in  agreement  with  predicted  ones 
based  on  observed  quantities  at  the  initial  point  where  a  path  was  straight. 
The  deviation  of  the  actual  paths  from  theoretical  ones  after  one  half 
wave  length  was  explained  as  an  effect  of  vorticity  transport  due  to 
lateral  mixing  process  (Ichiye,  1957).  Defant  (1956)  also  noticed  simi¬ 
larity  between  the  inertia  motions  in  the  equatorial  zone  studied  by  Whipple 
(1917)  and  Rossby's  meandering  trajectories  of  a  planetary  jet.  He 
speculated  that  the  eastward  flow,  such  as  the  Gulf  Stream  extension, 
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may  meander  with  large  amplitudes,  occasionally  discharging  eddies  and 
yet  remain  narrow  owing  to  its  stability,  as  shown  by  Rossby1  s  trajectories, 
while  the  westward  flow  like  the  North  Equatorial  Current  is  essentially 
unstable  as  eastward  trajectories  indicate  and  thus,  the  disturbances  are 
widely  spread  so  as  to  make  the  current  broad  by  lateral  transfer  of 
momentum.  Saint  Guily  (1957)  modified  R.ossby's  original  theory  by  including 
effects  of  variable  velocity  and  volume  transport  along  the  path  to  explain 
schematically  excessive  meanders  of  the  Gulf  Stream  as  observed  during 
Operation  Cabot. 

The  equation  of  a  trajectory  of  constant  absolute  vorticity  is 
derived  for  a  two-layer  system,  as  follows.  It  is  assumed  that  only  the 
upper  layer  with  a  depth  D  is  in  motion  and  the  lower  layer  is  infinitely 
deep.  Then,  the  potential  absolute  vorticity  is  constant  along  a  streamline 
(Charney,  1955)  or 

(£,  +  f)  D =  const.  (  on  a  streamline)  (1) 

where  4  is  the  relative  vorticity.  The  volume  transport  is  constant  along  a 
strip  in  the  upper  layer  bounded  by  two  neighboring  streamlines  and  the  inter¬ 
face.  Therefore,  we  have 

q  B  D  =  const  (2) 

where  q  is  the  velocity  and  E  is  the  width  between  the  streamlines  and 
a  bar  indicates  the  average  over  a  cross-section  of  the  strip. 

The  relative  vorticity  consists  of  curvature  and  shearing  terms, 


like 


^  =  lK  +  9q/9  n 


(3) 


where  K  is  the  streaxnline  curvature  and  n  is  a  coordinate  normal  to  the 
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streamline.  v'hen  two  streamlines  of  the  strip  are  close  enough,  the 
lateral  change  of  q  and  k  can  be  neglected  and  the  average  of  £  a 
over  a  cross-section  of  the  strip  or  jet  can  be  expressed  by 

t,  *  qK  ~  ~q  *~K  (4) 


In  the  following,  a  bar  denoting  the  average  is  dropped.  Then  after 
being  averaged  over  a  cross-section,  Equation  (1)  becomes 

D  (  f  +  q  •  K)  =  const,  (on  a  streamline)  (5) 

If  the  origin  of  x-y  coordinates  is  taken  at  a  point  of  a  straight  portion  of 
the  flow  and  x-  and  y-  axes  are  taken  in  the  east  and  north  direction 
respectively,  Equation  (5)  yields 


P  y  +  K  (  V/  D  E)  =  f0  (  D  /  D0  -  1 )  (6) 

where  the  suffix  o  denotes  the  value  at  y  =  0  and  V  is  equal  to  (DBq)Q 
or  the  volume  transport  of  a  jet  and  constant  along  the  jet  from  Equation 
(2).  Saint  Guily  (1957)  obtained  Equation  (6)  from  a  normal  component 
of  equations  of  motion  in  a  natural  coordinate  system, 

K  q2  +  fq  =  -  9p/  8n  (7) 

If  Equation  (7)  is  integrated  between  two  sides  of  the  jet,  the  righthand 
side  gives  pressure  difference  between  both  sides  and  is  constant  along 
the  jet  because  the  sides  consist  of  isobars.  'When  Equation  (2)  is  used, 
this  integrated  equation  becomes  (5), 

When  y(x)  denotes  the  coordinate  of  the  path  of  the  jet, 


curvature  K(x)  can  be  expressed  by 

K(x)  =  y"  i  1  +  (y')2 


-3/2 


where  a  prime  indicates  differentiation  with  xe  If  the  velocity  q  is 
constant,  D  is  also  constant  along  the  jet  and  Equation  (6)  becomes 


(8) 


55- 


in  which 


my"  +  y 

m  =  q  /  P 


1  +  (y')2 


-3/2 


(9) 

(10) 


Equation  (9)  can  be  solved  analytically  with  the  initial  conditions  that 


y(o)  =  0,  y'  (o)  =  tan  0 
The  first  integration  of  (9)  yields 


(11  a,  b) 


Y'  = 


(cos  0  +  y^/  2m)  2-l 


1/2 


(12) 


It  is  at  once  seen  that  the  ordinate  has  an  extreme 

yo  =  2m(l-cos0)  ^2 

which  is  the  amplitude  when  the  path  is  oscillatory.  Introducing  a  new 


(13) 


variable  4>  defined  by 

k2  sin2  4-  =  (  y^  /  2  m  )  (  1  +  cos  4>  +  y2/  2  m) 
k2  =  (  1  -  cos  (  )  /  2 


-1 


with 


(14) 

(15) 


we  can  integrate  (12)  which  becomes 

x  =  n/  m  cos  0  ,  (  1  +  k2  sin2  <l/ cos  0)  (  1-k2  sin24>^2  d  4>  ( 1 


J 


o 


or 


x=^m!2E(9,k)-F(4>,  k)-k2  sin  (  1  -k2  sin^  )  ^  ^ 


(17) 


where  F  (4>,  k)  and  E  ($,  k)  is  Legendre's  elliptic  integral  of  tha  first 
and  second  kind,  respectively.  An  extreme  yQ  corresponds  to  4>  =  tt  / 2 


or  the  value  of  x  defined  by 

o  1 


x  =  Njm  (2  E  -  K) 


(18) 


where  E  and  K  are  complete  elliptic  integrals.  The  paths  for  0  =  45® 
90°,  117°  30',  180c  are  shown  in  Fig,  8,  Then  0  reaches  117®  30',  the 
path  touches  itself  (occlusion)  making  looplike  meanders,  V.’hen  0  exceeds 
130°  44',  the  jet  flows  westward,  making  a  series  of  loops  and  at  0  =  180°, 


the  path  forms  only  one  loop  at  x  =  -c>Q.  The  values  of  xq  (a  quarter  of  a 
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wave  length)  yQ  (half  amplitude)  for  different  values  of  6 


6  (r) 

15 

45 

75 

90 

117c  5 

o/N'm 

0,26 

0.76 

1*22 

1  e  4  1 

1,71 

1. 55 

1.39 

1  o  07 

0C  85 

0,  32 

alues  of 

m  in  km. 

are 

for  5,  25, 

and  100 

cm/ sec  of 

respectively,  47,  105  and  210  at  the  Equator  and  55,  120  and  240  at  45' 
latitude,. 

In  order  to  discuss  the  effect  of  change  of  depth  D  on  the 

with 

course  of  meanders,  it  is  assumed  first  that  q  varie\the  distance  2 
along  the  jet  from  an  origin  for  constant  D  and  B„  Then,  Equation  (6) 


oecome  s 


K  q  (  2  )  +  (5  y  =  0  (19) 

When  q  decreases  v/ith  2  ,  the  curvature  K  increases,  the  angle  of  the 
path  with  y»  axis  0  increases  and  steepness  (ratioiof  amplitude  to  wave 
length)  of  meanders  increases,  leading  occasionally  to  occlusion,,  On  the 
other  hand,  when  q  increases  with  2  ,  K  decreases,  0  decreases,  the 
steepness  decreases  and  the  path  of  the  jet  is  more  straight,  Nexts  it  is 
assumed  that  the  depth  D  is  dependent  on  x  only  and  that  B  is  constant. 
Then,  if  y  is  replaced  by  yj  =  y  ♦*  ( fQ /p )  (  D/  Dq  »  1)  in  Equation  (6),  the 
same  conclusion  as  above  will  be  obtained,  i0e0  steepness  of  meanders 
increases  when  the  depth  increases  with  x  and  vice  versa.  The  only  dif¬ 
ference  is  that  the  axis  of  the  meanders  (the  mean  path)  is  y,  =  0  in  this 
case  instead  of  y  =  0  (x«"  axis)  for  the  case  of  constant  q,  When  the 
depth  is  a  linear  function  of  y  like  D  =  D0  (  1  +  sy),  Equation  (10) 
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yieids 

Kqo  (  1  +  sy)'1  +  (p  -  sfQ)  y  =  0  (20) 

when  s  •=  sQ  =  pf0“‘,  the  p-effect  is  compensated  by  the  effect  of  depth 

and  the  jet  flows  straight,,  as  noticed  by  Ekrnan  (1923)  and  Neumann  (1955), 

From  Equation  (20),  a  general  rule  can  be  obtained  according  to  s  <  sq 

or  s  >  s  £  '’’’hen  s  <  s  ,  steepness  of  meanders  is  less  for  a  >  0  and 

o  ° 

more  for  s  <  0  than  for  a  constant  depth,  When  s  >  s  and  s  >  0,  an 

o 

eastward  jet  becomes  unstable  and  an  occlusion  occurs  immediately,  while 
a  westward  jet  shows  stable  meanders.  When  |  s|>  sq  and  s  <  0,  vice 
versa, 

Ichiye  (1957)  extended  Equation  (6)  to  a  case  that  the  initial 
angle  9  and  the  velocity  q  vary  periodically  with  time  and  explained  the 
coupling  of  a  cyclonic  detour  to  the  south  of  the  Kuroshio  between  135®  and 
140°  E  and  its  northward  anti-cyclonic  meander  east  of  142°  E  from  an 
equation  modified  from  (6),  Also,  he  considered  the  horizontal  transport  of 
relative  vorticity  across  both  sides  of  a  jet  through  eddy  viscosity,  approxi¬ 
mated  it  with  a  linear  function  of  a  distance  SL  and  equated  it  with  the  right- 
hand  side  of  (6).  It  is  shown  that  the  influx  of  positive  vorticity  causes  less 
cyclonic  curvature  and  farther  northward  extension  of  the  current,  suggest¬ 
ing  that  an  occasional  occlusion  of  the  Kuroshio  east  of  150c  E  is  due  to 
lateral  influx  of  relative  vorticity, 

Warren  (1962)  studied  the  effects  of  bottom  topography  on  the 
course  of  the  Gulf  Stream  north  of  Cape  Hatteras,  From  conservation  of 
absolute  vorticity,  he  derived  the  equation  which  is  essentially  the  same 
as  Equation  (6)  excepting  that  D  is  equal  to  the  actual  depth  instead  of  the 
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depth  of  the  upper  layer.  He  considered  an  example  that  D  is  a  linear 
function  of  both  x  and  y  like 


D-D^+SjX+s^y 

Further,  instead  of  a  simple  form  of  velocity  q  and  volume  transport 


(21) 


q  D  B,  he  used  quantities  integrated  across  the  entire  stream  width,  such 


as 


C  2, 


V  =  \  q  D  d  n,  M.  =  \q  Ddn,  T  =  \  q  d  n 


^  J*  '  J 

Then,  the  averaged  vorticity  equation  similar  to  (9)  becomes 


where 


(22) 


K  (x)  +  Njy  -  N2x  =0  (23) 

Nx  =  (PV  -  fQTs2)  ,/T1,  N2  =  fQ  Tsj  M**1  (as  a,  b) 

By  rotation  of  x-y  and  transformation  of  (x-y)  system  into  (x*  ~y')  system, 


Equation  (23)  is  changed  into  the  same  form  as  Equation  (9)  except  that 
m~2  =  N  2  +  No2 


(24) 


instead  of  (9).  The  analytical  solution  of  the  transformed  equation  is  given 
by  relations  (16)  or  (17)  with  x  and  y  replaced  by  x'  and  y'.  In  order 
to  determine  the  path  of  the  Gulf  Stream  east  of  Cape  Hatteras,  warren 
(1963)  applied  this  equation  to  a  series  of  small  segments  starting  at  a 
point  of  no  curvature.  He  used  the  observed  value  of  velocity  and  direction 
of  the  Gulf  Stream  at  a  starting  point  only  and  successively  proceeded  to 
each  segment  in  which  s^  and  s2  are  determined  from  a  chart.  One  of 
his  results  which  showed  the  best  agreement  with  an  observed  path  is  shown 
in  Fig,  8,  On  some  occasions,  observed  paths  showed  larger  meanders, 
than  calculated,  suggesting  that  an  eddy  influx  of  vorticity  across  the  sic ..  j 
may  be  important  as  in  the  case  of  the  Kuroshio,  discussed  by  Ichiye  (1957). 
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Greenspan  (1963)  also  noticed  that  the  10-,  100-  and  1000- 
fathom  lines  almost  coelesce  just  south  of  Cape  Hatteras  and  suggested 
that  the  Gulf  Stream  being  stable  as  an  inertial  boundary  current  south 
of  the  Cape  becomes  unstable  owing  to  a  sudden  increase  of  the  depth 
north  of  it.  His  argument  was  based  on  his  criterion  for  the  existence 
of  inertial  boundary  layers  (Greenspan,  1962).  In  a  homogeneous  ocean 
with  a  variable  depth  D  or  in  a  two  layer  system  with  an  upper  layer 
of  a  thickness  L  and  with  a  motionless  lower  layer,  conservation  of 
volume  and  absolute  vorticity  yields  two  equations  (Charney,  1955). 


uD  =  -9vJj/9y,  vD  =  9i|j/9x 

(25) 

(f +U/D  =  F  U) 

(26) 

and 


where  F  is  an  arbitrary  funetion\  t,  is  the  relative  vorticity.  As  in 
Section  6,  these  equations  can  be  expressed  with  non-dimensional  vari¬ 
ables  by  introducing  U,  L,  H  and  f  as  characteristic  velocity,  length 
depth  and  Coriolis'  coefficient,  respectively.  Then,  the  absolute  vorticity 
is  given  by 

f  +  €  V  *  (  D"*1  V  ^  )  (27) 

where 

<=  =  U/(2f0L)  (28) 

is  a  small  dimensionless  parameter  if  L  is  taken  as  the  lateral  dimension 
of  the  ocean.  In  Equation  (27)  and  the  following,  all  the  quantities  are 
expressed  in  a  dimensionless  form.  A  new  coordinate  system  £=£  (x,  y), 
q  (x,  y)  is  introduced  by  use  of  a  conformal  transforxnation 

4  +  irj  =  w  (x  +  i  y )  (29) 

with  w  an  analytic  function.  Furthermore,  let  £  =  £  B*  ^  <  and  -q 

denote  the  coast,  the  ocean  and  a  coordinate  normal  to  the  coast,  respectively. 


, 


—  GO i* 


The  stream  function  ijj  varies  rapidly  across  the  narrow 
width  of  a  coastal  boundary  layer  from  its  interior  value  T  (£>r)  )  to 
zero  at  the  coaste  Therefore?  by  introducing  a  stretched  boundary 
coordinate 

Z  =  (6  -  £B>  6  ’a  (  a  >  o  )  (30) 

The  stream  function  may  be  divided  into  two  parts: 

4*  =  $  (z»r|)  +  T|)  (31) 

Transforming  Equation  (26))  substituted  by  (27)  by  use  of  (29)  and 
(3ofj  and  equating  the  coefficients  of  e  “1/2  £or  a  =  1  / Z  to  zero,  we 
have 

4>  z2(JD)“2  +  fD“1=C(^H)  (32) 

where 

J2  =  x  ^  +  y^Z  (33) 


and  Q  is  an  arbitrary  function.  If  an  inertial  boundary  layer  is  to  exist. 


then  4>  must  vanish  at  the  seaside  edge  of  the  boundary  layer  (z  =  -  cq  )  ♦ 


Therefore, 


(  T  )  =  f  D 


-1 


(34) 


with  T  =  \3>  (  £  ,  r)  ).  Equation  (32)  becomes 

jD 


♦  ZZ  =  (JE)' 


0(4+  <£  )  -  0  (  4 ) 


"I 

J 


(25) 


For  large  z,  4>  is  so  small  that  the  bracketed  term  can  be  expanded  as 


a  Taylor  series.  Thus, 

$  *  (  J  D)2  (  9  Q  /  9  ^  )  <t> 


zz 


(36) 


v/hen  z  ~  -  oo  •  However,  <J>  approaches  asymptotically  zero  as  z 
increases  only  if  the  coefficient  of  <j>  in  (36)  is  positive  (  $  becomes 
oscillatory  with  z  if  the  coefficient  is  negative)  .  Therefore,  the  condi¬ 
tion  for  existence  of  the  boundary  layer  is  simply 


ac  /  9^  >  o 


(37) 


4  1 


C-» 


for  |  ~  By  use  of  (24),  tbj.3  condition  becomes 

9  (  f  /  D)  ^  >0  (38) 

However,,  ^  is  the  stream  function  of  the  interior  flow  and  thus  - 

■n 

is  tbo  |  component  of  the  interior  flow  in  the  direction  of  increasing 


—4 

rj  »  If  Vj  and  n  denote  the  vectors  of  the  interior  flow  velocity  and 
the  unit  vector  normal  to  the  boundary  in  the  outward  direction,  we 
have 


— t  — > 

Si  =  -  vT  •  n. 
r\  I 

Therefore,  condition  (38)  can  be  written  as 

y  _ _ y 

(  Vj  •  n  )  a(f/D)/Srj<0  at  £  =  |B 


(39) 


Greespan  (1963)  concluded  that  the  effect  of  topography  south 
of  Cape  Hatteras  controls  the  path  of  the  Gulf  Stream  as  an  inertial  boundary 
current  since  the  criterion  (39)  is  satisfied  there  owing  to  a  slight  change 

—  y  — 

of  D,  while  9  f  /  9  p  >0  and  Vj  «  n  <  0.  Just  south  of  the  Cape,  the 
positive  value  of  9  D/  9t|  exceeds  the  effect  of  the  northward  increase 
of  f  and  the  criterion  (39)  is  violated,  He  suggested  that  *he  control 
of  topography  may  be  effective  even  to  a  two  "layer  system  through  some 
intermediary  means,,  Also,  he  speculated  that  the  location  and  flow 
structure  at  the  separation  point  just  couth  of  the  Cape  may  change  with 
the  order  cf  several  weeks  according  to  criterion  (39),  since  in 
■Warren's  calculation  (1962)  the  initial  condition  for  the  Gulf  Stream 
somewhere  just  north  of  the  Cape  must  be  taken  differently  for  each 
meander. 
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80  Energy  Transfer  Between  the  vleanders  and  the  Mean  Flow 

Since  vorticity  transfer  through  the  sides  of  the  flow  seems  to 
play  an  important  role  in  meanders,  as  discussed  in  Section  7,  Ichiye^  1957) 
tried  to  determine  momentum  transfer  in  the  Kuroshio  due  to  meanders.  The 
G  E  K  data  obtained  in  the  area  between  135°  E  and  150°  E  in  March  and 
May,  1956  were  divided  in  two  groups  east  and  west  of  longitude  141°  E, 
each  group  containing  data  on  seven  to  nine  transects  across  the  currents 
consisting  of  stations  of  less  than  ten  miles  apart  along  each  transect.  In 
each  transect,  the  current  at  each  station  was  decomposed  into  components 
parallel  and  normal  to  the  mean  flow  over  the  transect.  Then  each  component 
v/as  grouped  again  according  to  the  distances  from  the  axis  (having  maximum 
speed)  of  the  current,  and  averaged.  Also,  for  each  group,  correlations  of 
deviations  of  velocity  u*  and  v*  were  computed,  where  u  and  v  are 
velocity  components  normal  and  parallel  to  the  mean  flow  determined  for  each 
transect.  Therefore,  if  a  bar  denotes  the  average  over  a  group  consisting  of 
stations  of  equal  distance  from  the  current  axis,  the  energy  transfer  from 
meanders  or  (eddies)  to  the  mean  flow  is  given  by 

p  u'  v'  37/  9x  (1) 

where  x  denotes  the  distance  normal  to  the  mean  flow.  The  computed  values 
of  this  transfer  and  mean  speed  u  for  the  western  and  eastern  region  are 
shown  in  Fig.  10-a  .  In  the  western  region,  the  energy  is  transferred  from 
the  mean  flow  to  meanders  to  the  left  side  of  the  axis  (facing  downstream) 
and  in  the  reverse  direction  to  the  right  side  but  in  the  eastern  region  vice 


versa. 
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v  ebster  (1961a)  computed  the  energy  transfer  from  eddies 

to  the  mean  flow  at  two  transects  across  the  Gulf  Stream  from  GEK  data, 

using  a  relation  similar  to  (1)  except  the  bar  denoting  a  time  average  of  a 

quantity*  One  transect  is  between  Tiami  and  Gun  Cay,  Bimini.  A  total 

of  632  GEK  readings  from  42  cruises  made  from  1956  to  1958  were  divided 

into  eleven  longitudinal  zones  and  the  average  shear  9  v/  8  x  and  Reynolds 

stress  u'v'  are  computed  for  each  zone.  Another  transect  off  Onslow 

Bay  was  crossed  almost  normal  to  the  Gulf  Stream  axis  one  hundred  and 

twenty  times  during  some  thirty  consecutive  days.  The  GEK  readings  were 

divided  into  eleven  zones  and  treated  in  the  same  ways  as  Florida  Straits 

data.  The  results  shown  in  Fig.  10b  indicate  that  the  energy  is  transferred 

from  meanders  to  the  mean  flow,  contrary  to  the  conventional  speculation 

that  the  meanders  represent  a  mechanism  for  frictional  dissipation  of  the 

mean  flow.  This  feature  is  similar  to  the  energy  transfer  from  large  scale 

eddies  to  the  mean  zonal  flow  in  the  atmosphere  determined  by  Starr  (1953), 

whose  results  are  shown  in  Fig,  10c  by  plotting  the  production  of  zonal  kinetic 

energy  and  the  relative  angular  velocity  as  a  function  of  latitudes.  Though 

the  energy  transfer  has  a  similar  feature  between  two  transects,  its  values 

averaged  over  the  entire  length  of  each  transect  are  quite  different  and 

equal  to  79  x  10“^  and  3  x  10~^  ergs  cm  3/  sec  for  off  Onslow  Bay  and 

the  Florida  Straits,  respectively.  Since  the  kinetic  energy  of  the  mean 

3  3 

surface  flow  averaged  across  the  stream  is  7.3  x  10  ergs/ c.m'5  and 
•2 

8C  6  x  10  for  off  Onslow  Bay  and  the  Florida  Straits,  respectively,  the 
calculated  rate  of  energy  transfer  doubled  the  mean  surface  kinetic  energy 
in  1 1  in  the  former  and  329  days  in  the  latter.  Webster  also  discussed 


. 
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the  energy  equation  for  the  mean  flow  in  a  steady  state 

^  ^ 
7h 


~ >  —  =T  -  __ 

V  •  Vjj  E  +  u  V  •  p  u'v'  +  v  V  *  p  v'v'  +  V*VVip+V*  F  =  0 


which  can  be  obtained  from  equations  of  motion  in  a  horizontal  plane  and 

— b’  —  —  — y 

V  =  ( u,  v),  v1  =  (u'»  v1,  w'),  ^  =  ( 3  /  9  x,  3/3  yi  ?  V  =  (  3/3x,  3/3  y, 


9  /  3  z)  ^  E  =  (  u^  +  v^)/  2  and  F  is  the  frictional  forcee  Integrating  (2) 
over  a  cross-section  of  the  current,  he  estimated  that  the  contribution  to 
the  kinetic  energy  of  the  mean  flow  by  horizontal  motion  represented  by 
the  sum  of  the  second  and  third  term  of  (2)  comes  mainly  from  the  term 
given  by  (1)  while  the  effect  of  wind  stress  expressed  by  the  last  term  of 
(2)  is  less  than  the  contribution  from  (1)  by  one  to  two  orders  of  magnitudes. 

Webster  (1961b)  analyzed  the  data  of  the  transects  off  Onslow 
Bay  and  found  that  the  amplitudes  of  harmonic  analysis  of  the  location  of  the 
20°  isotherm  at  200  m  depth  and  of  the  velocity  maximum  determined  by 
GEK  along  the  transects  have  peaks  at  4  and  7  days  with  values  of  6  miles 
and  3,7  miles,  respectively.  He  concluded  that  there  is  no  definite  evidence 
for  meanders  off  Onslow  Bay  as  relating  to  the  tidal  period,  suggesting  that 
they  may  be  related  to  local  meteorological  disturbances.  He  proposed  a 
model  of  meanders  in  this  area  consisting  of  those  associated  with  an 
intense  offshore  current  with  one  to  four  days'  duration  and  those  accompany* 
ing  a  broad  onshore  current  of  four  to  seven  days'  duration  and  of  a  wave 
length  of  a  hundred  kilometers  in  contrast  with  long  meanders  north  of  Cape 
Hatteras  having  wave  lengths  of  1000  km  or  more  and  periods  of  several 


months. 


Oort  (1963)  computed  co-variance  in  time  of  simultaneous  cross 


(2) 


current  velocity  (  u' )  temperature  (  T' )  and  density  (  p *  ) ,  u'T'  and  u'p' 
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from  three  transects  including  the  one  off  Onslow  Bay  analyzed  by  V  ebster 
(1961a).  Other  transects  are  70  crossings  along  the  30°  N  parallel  off 
Jacksonville,  in  October -November,  1961  and  20  crossings  along  lines 
shifting  slowly  downstream  off  Cape  Hatteras  in  October,  1962.  He  used 
surface  data  of  temperature,  density  and  GEK  readings,  as  well  as  BT 
data  for  temperature  at  30  m,  60  m  and  100  m  with  surface  GEK  data 
for  currents.  It  was  found  that  the  effect  of  a  rotation  of  the  coordinate 
axes  by  about  10°  did  not  change  the  order  of  magnitude  and  sign  of  u'T' 
suggesting  that  v'T'  is  amall  compared  with  u'T'.  Also,  since  a  slight 
rotation  of  axis  may  make  the  mean  cross  current  velocity  u  vanish,  the 
total  heat  transport  across  the  stream  uT  can  be  well  approximated  by 
the  eddy  heat  transport  u'T'.  The  heat  transport  (u'T1)  ^temperature 
and  the  downstream  velocity  at  the  surface  at  three  transects  are  shown 
in  Fig.  11.  The  results  at  three  transects  show  that  the  offshore  heat  flux 
coincides  approximately  with  the  place  of  maximum  current.  This  may 
explain  the  appearance  of  a  warm  core  to  the  right  of  the  current  peak 
(facing  downstream),  which  is  apparent  in  the  data  shown  in  this  figure. 
Also,  these  data  indicate  heat  fluxes  from  the  colder  coastal  waters  to 
the  warmer  ocean  in  all  three  transects.  The  surface  heat  flux  aver¬ 
aged  along  three  transects  is  1 . 8  ±  1 . 3,  3,  5  ±  1 . 6  and  0,46  ±  0,55 
(°C  cm/ sec),  respectively  off  Cape  Hatteras,  Onslow  Bay  and  Jacksonville. 

The  average  density  flux  u'p'  is  -0.9  ±0.5,  -1,8  ±0.6  and  -0,15  ± 

-3-1-2  ... 

1,9  (10  g  sec  cm  )  from  north  to  south,  again  indicating  the 

mass  transport  against  the  horizontal  density  gradient.  The  heat  flux 
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at  deeper  layers  shows  the  same  tendency  as  the  surface  values  and  its 
magnitude  averaged  on  a  transect  is  three  to  ten  times  the  surface  value 
off  Jacksonville*  However,  values  of  the  subsurface  heat  flux  is  subject 
to  more  errors  due  to  the  use  of  surface  currents  instead  of  currents  at 
corresponding  depths^  since  a  strong  vertical  shear  usually  exists  near 
the  surface.  An  idea  of  the  order  of  magnitude  of  computed  heat  flux  may 
be  gained  by  calculating  the  eddy  heat  flux  normal  to  a  vertical  wall  along 
the  Gulf  Stream  for  a  distance  of  1000  km  and  with  a  depth  of  100  m  the 

heat  transport  p  Cu1  T1  across  the  wall  (  C,  heat  capacity  of  sea  water) 

18  l  - 1 

becomes  10  cal  day”1  corresponding  to  u'T*  =  10  (°C  cm  sec  ),  This 

value  is  comparable  with  the  calculated  annual  poleward  heat  transport 
18  —  1 

of  20  x  10  cal  day  at  30°  N  along  a  distance  of  6700  km  of  the  entire 
North  Atlantic  Ocean,  Therefore,  the  meanders  play  an  important  role 
in  maintenance  of  the  Gulf  Stream  through  thermal  process  as  well  as 
through  momentum  transport  discussed  by  Webster, 

Chew  (1958)  proposed  a  hypothesis  that  meanders  due  to 
internal  waves  of  a  semi-diurnal  tidal  period  may  cause  fluctuations  of 
cross  stream  components  of  surface  currents  measured  with  GEK  between 
Miami  and  Cat  Cay,  Chew  (1961)  also  determined  the  average  surface 
cross  stream  component  on  this  section  from  500  GEK  data  over  eight 
years  and  1100  data  from  Pillsbury' s  measurements.  The  distribution 
indicates  a  westward  (offshore)  flow  for  a  distance  of  20  km  from  the 
Florida  Coast  with  maximum  speed  of  10  c,m /  sec  and  an  eastward 
(onshore)  flow  to  the  east  until  Cat  Cay  with  the  maximum  speed  of 


' 
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14  cm/ sec  at  about  79a40,E.  He  attributed  the  divergence  near  the  coast 
to  the  lateral  Ekman  layer  due  to  the  coast  and  lateral  eddy  viscosity.  It 
is  yet  to  be  seen  how  this  systematic  distribution  of  the  mean  cross  current 
component  may  influence  the  momentum  and  heat  transfer  across  the  stream 
determined  by  Webster  (1961a)  and  Oort  (1963). 
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9.,  Model  Experiments  for  Meanders 


A  model  experiment  of  the  wind-driven  ocean  circulation  by 
use  of  a  rotating  circular  tank  was  initiated  by  Von  Arx  (1952),  who 
constructed  a  model  of  the  entire  northern  hemisphere  in  a  paraboloidal 
basin  of  two  meters  in  diameter  and  of  a  focal  length  of  0,5  meter.  His 
results  showed  a  strong  similarity  to  the  actual  circulation  both  in  the 
North  Atlantic  and  North  Pacific  Ocean  to  details  of  the  meanders  of  the 
Gulf  Stream  between  Cape  Hatteras  and  the  Grand  Banks  and  of  the 
cyclonic  bend  of  the  ICuroshio  off  the  southern  coast  of  Japan,  The  occur¬ 
rence  of  a  western  boundary  current  like  the  Gulf  Stream  and  its  meander¬ 
ing  is  mainly  due  to  the  (5-effect,  as  discussed  by  many  authors  and  it 
was  found  that  the  similarity  in  the  (5-effect  between  the  prototype  and 
the  model  may  be  established  by  changing  the  depth  D  of  a  model 


according  to  the  relationship 


(8  £/  9y)  =  fa  U  8y  +  (  f  +  £,)  D-1  8  D/  8  y 

P  I 


s 


m 


(1) 


where  subscripts  p  and  m  refer  to  prototype  and  model,  respectively; 
y  is  the  poleward  coordinate  (or  the  radial  coordinate  in  the  model),  t, 
is  the  relative  vorticity.  Further,  since  for  a  flat  basin  (8  f  /  9  y)  =  0 
and  also,  £,  is  usually  negligible  compared  with  f,  Equation  (1)  indi¬ 
cates  that  the  constant  increase  of  depth  D  of  the  model  towards  the 
rotation  center  according  to  (1)  may  simulate  the  (3-effect,  A  rotating 
model  experiment  with  a  1, 5  meter  diameter  basin  made  by  Ichiye  and 
Plutchak  (19  62)for  the  circulation  in  the  Gulf  of  Mexico  showed  that  the 


modification  of  depth  according  to  (1)  produced  an  anti-cyclonic  loop  of 
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the  flow  after  entering  the  Gulf  through  the  Yucatan  Straits*  as  observed 
in  hydrographic  surveys. 

The  occurrence  of  meanders  in  the  eastward  jet  was  also 

noticed  by  a  qualitative  experiment  in  a  pie-shaped  basin,  where  an  ocean 

circulation  similar  to  the  real  one  was  produced  by  suitable  distribution 

of  a  sink  and  a  source  (Stommel  and  others,  1958),  Faller  (1961)  found 

that  the  compensation  flow  in  the  bottom  Ekman  layer  combined  with  a 

radial  depth  variation  in  the  model  yields  an  equation  of  a  stream  function 

similar  to  the  equation  for  a  wind-driven  ocean  circulation,  A  slight 

modification  of  his  analysis  is  summarized  as  follows: 

For  a  fluid  in  a  rotating  model  basin,  r,  0  and  z  are  taken 

radially  outward^in  the  direction  of  rotation  (counter-clockwise)  and 

vertically  upward,  respectively  and  us  v  and  w  are  velocity  components 

respectively  in  r~,  0-  and  z-  directions,  When  the  Rossby  nuxnber  of 

the  system  is  sufficiently  small,  the  horizontal  velocity  components  can 

be  represented  by  geostropnic  relations 

u  =  u.  =  -  r1  9h/(  rae),  v  =  v  =  f"  !ah/ar  (2  a,  b) 

8  § 

where  h  is  the  depth  outside  the  bottom  Ekman  layer  with  depth  De  = 

v  n/  (  2  v!  f)  in  which  1/  is  a  molecular  viscosity.  At  the  bottom  z  =  0, 

velocity  and  its  vertical  shear  should  vanish.  Since  De  <<  h,  the  condi- 

at 

tion  of  a  geostrophic  flow \  z  -  h  yields  the  horizontal  velocity  as 

u  =  ug  -  exp  (-x)  ( ug  cos  x  +  vg  sin  x)  (3) 

v  —  Vg  —  exp  (-x)  (Ug  sin  X  “  vg  cos  x)  (4) 

where  x  =  nz  /  De.  The  equation  of  continuity  with  a  source  P  per  unit 


of  volume  is 


..  . 
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r“  *  9(ru)/  9  r  +  9  v  /  (  r  9  0)  +  9w/9z  =  P 


(5) 


Substituting  (3)  and  (4)  into  (5),  we  have 

8  w/3  z  +  exp  (  -  x  )  sin  ^  j^3  u^  /  3  0  -  9(  r  v  )!  9  r  ]  r"^=  P  (6) 
In  the  following,  subscript  g  is  dropped.  At  the  surface  z  =  h,  the 
boundary  condition  for  w^  becomes 

w^  =  u  (3  h/  3  r)  +  v  3  h/  (  r  3  6)  +  9  h/  3  t  (7) 

where  3  h/  3  t  is  change  of  depth  due  to  steady  vertical  motion  of  the  free 
surface. 

Integrating  (6)  with  z  from  o  to  h  and  using  (7),  we  have 


u3h/ 3  r  +  v  3  h  /  (  r  9  0}  +  3h/3t  +  De 
-  3(rv)/3r]  /  (2  tt  r)  =  S 


3  u/  9  e  - 


J 


where  S  =  i  P  d  z, 
^o 


Introducing  a  stream  function  defined  by 


(8) 


u  =  3  4*  /  (  r  3  ©),  v=~34>/0r  into  Equation  (8),  we  have 


(3  h  /  3  r)  (Svjj/  r  3  0}  +  De  (Z  tt)  ” 1  V2^  =  S  -  3h/  9t 


(9) 


This  equation  is  similar  to  the  equation  for  a  wind-driven  circulation 
derived  by  Stommel  (1948)  if  3  h/  3  r  and  S  -  3  h/  3  t  are  considered 
as  p  and  the  wind  stress,  respectively. 

Simulation  of  the  p-effect  by  changing  a  depth  may  be  achieved 
either  by  use  of  a  distorted  depth  distribution  or  by  taking  advantage  of  rise 
of  the  free  surface  toward  the  rim  due  to  a  centrifugal  effect.  The  former 
approach  was  taken  by  Ichiye  and  Plutchak  (1962)  for  the  Gulf  of  :dexico 
experiments.  The  latter  approach  was  made  by  Stommel  and  others  (195  8) 
for  their  qualitative  experiments  in  a  pie-shaped  basin.  The  equation  of  the 


. 
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free  surface  is  given  by 

h  =  h0  (  1  +  £2  r2/4  gho  )  (10) 

where  a  is  the  radius  of  the  basin  and  hQ  is  the  depth  at  the  center  of  the 
basin.  Therefore,  by  taking  suitable  values  of  rotation  speed  and  initial 
depth,  the  ^--effect  may  be  simulated  for  a  flat  bottom  tank. 
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(ADDENDA) 

A  series  of  experiments  including  determination  of  stability  of 
various  kinds  of  jets  and  generation  of  vortices  due  to  meanders  are  now 
underway  by  N,  B.  Plutchak  and  the  author,  using  a  rotating  basin  of  1.5 


m  diameter. 
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Explanation  of  Figures 


Fig.  1 

Positions  of  the  Gulf  Stream  at  the  beginning  and  end  of 
Operation  Cabot  after  Fuglister  and  "Worthington  (1951). 
Isotherms  of  mean  temperature  (®F)  of  the  upper  200- 
meter  layer  and  belts  of  the  "warm  core"  are  shown 
west  and  east  of  64°  W,  respectively. 

Fig,  2 

(A)  Temperature  (°C)  at  the  200-meter  layer  from 

Gulf  Stream* 60  (May-June,  I960) 

(B)  Positions  of  the  abrupt  change  in  direction  of  the 

Gulf  Stream  and  the  position  of  the  maximum  gradient 
in  the  average  200-meter  temperature  (after  Fuglister 
(1963)}  L 

Fige  3 

Velocity  distributions  in  Haurwitz  and  Panof sky’s 
model  (a,  b  and  c)  and  in  Kuo's  model  (d).  (Numerical 
values  in  (d)  are  those  used  for  Ichiye's  calculation 
of  Kr  and  Ki  in  Fig.  4). 

Fig.  4 

Roots  of  the  frequency  equations, 

I  and  II  (  H-P  theory  ):  CH  =  (  c  -  Uj)  /  (  U2  3  -  UL  ), 

d2  =  100,  d3  =  20  d  »L,  d2  =  20  (for  I)  and  200  (for  II) 
all  in  nautical  miles. 

K  (Kuo’s  theory):  CK  =  (c-U0)  /  Um  -  UQ),  Um  =  1.  1  (kt) 
UQ  =  -  0a  1  (kt)  (See  Fig.  3~d) 

Suffices  r  and  i  denote  respectively  real  and  imaginary 
part  . 

Fig.  5 

Diagram  of  the  two-layer  model  used  by  Iida  and  Ichiye  . 

Fig.  6(a) 

Diagram  of  the  two-layer  ocean  used  by  Duxbury:  the 
equilibrium  state  is  indicated  by  full  lines,  the  disturbed 
state  by  dashed  lines.  The  plane  y  =  0  represents  a 
velocity  shear  boundary  in  the  equilibrium  siate 

(6b) 

Graph  of  the  negative  real  values  of  the  relative  frequency 
V  versus  the  relative  wave  number  p.  for  F"  =  10„0, 

The  mode  number  represents  the  numbers  of  zeros  of 
the  function  Z(y)  in  the  Region  1.  (After  Duxbury,  1963). 

(6c) 

Iwata’s  model  of  the  two-layer  system.  Both  sides  of  the 
uniform  jet  U  are  motionless  and  have  constant  depths 

L>2  and  D^. 

Fig.  7 

Semi»circle  in  complex  velocity  plane  in  which  the  phase 
speed  for  unstable  baroclinic  waves  must  lie. 

Explanation  of  Figures  ( coni' d) 


Fig.  8 

Trajectories  of  the  constant  absolute  vorticity  flow, 
after  Rossby  (1940)  and  Saint-Guily  (1957). 

(a)  I  :  0Q  =  45°,  II  :  6  =  90°  ,  III  :  eQ  =  1 17.  5° 

(b)  0  =  180° 

Fig.  9 

Gulf  Stream  path  observed  on  June  8,  1950  with  a  path 
calculated  by  applying  Equation  (23)  to  small  segments 
A  to  L.  After  Warren  (1963)  , 

Fig.  10 

Kinetic  energy  transfer  from  eddies  to  mean  flow 
(solid  line)  and  mean  downstream  velocity  (dotted  line) 

(a)  The  Kuroshio,  after  Ichiye  (1957),  (a-1):  Between 
135°  E  and  141°  E,  (a~2):  between  141°E  and  150°  E  . 

(b)  The  Gulf  Stream  after  Webster  (1961). (b-1):  Off 
Onslow  Bay,  (b~2):  Miami-CatCay 

(c)  The  upper  atmosphere  after  Starr  (1953).  (Dotted 
linerrelative  angular  velocity.) 

Fig.  11 

The  cross  current  eddy  heat  flux  (u'T1),  the  mean 
temperature  (T)  and  downstream  velocity  (v)  as  a 
function  of  distance  from  the  mean  axis  of  the  Gulf 
Stream.  A  positive  value  of  eddy  heat  flux  repre¬ 
sents  an  offshore  transport  of  heat 

(a)  Cruise  off  Onslow  Bay  from  May  20  to  June  18, 
1958. 

(b)  Cruise  off  Jacksonville  (30°N)  from  October  7  to 
November  15,  1961. 

(c)  Cruise  off  Cape  Hatteras  during  October,  1962 
after  Oort,  1963) 

Fig.  12 

Surface  dynamic  topography  referred  to  1000  d  bar 
level  of  the  Kuroshio  (a  and  b,  after  Ichiye,  1957a) 
and  bottom  topography  of  the  Kuroshio  area  (c,  after 
Ichiye,  I960) 

May-June,  1951  (Nc^feouthward  detour) 
b)  May,  1955  (Southward  detour  between  135°  and  1408  E) 
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Explanr-tion 

Fig.  13 


Fig.  14 


Fig.  15 


of  Figures  (cont'd,) 


(a)  Southern  limits  of  the  Kuroshio  between  135°  and 
140°  E  (Crosses  denote  the  values  determined  from 
dynamic  topography#  and  solid  circles  those  from 
isotherms  at  100-meter  layer)  and  appearance  periods 
of  an  isolated  cold  water  (hatched  bar)  after  Ichiye 
(1957). 

(b)  The  positions  of  northern  boundary  of  the  Kuroshio 
along  144°  E  (a  vertical  bar)  and  the  trend  curve  of 

coastal  temperature  varying  with  a  period  of  4e  5  years 

(after  Kawai ,  1955) 

The  positions  of  the  axis  of  maximum  current  in  the 
Kuroshio  determined  from  GEK  measurements  and 
subsurface  temperature  distributions  (a)  1956  to  1958 
(b)  1959  (after  Masuzawa,  1961) 

The  bathymetri  chart  (a)  and  temperature  and  current 
distributions  (b) 

(a)  based  on  Monaco  chart 

(b)  May-June,  1962.  Thick  arrows  indicate  surface 
currents  between  1.0  and  1.5  knots  and  thin  arrows 
between  0.5  and  0,9  knots  measured  with  GEK, 
Isotherms  at  100  meter  depth  are  in  centigrade, 

(after  Japanese  Hydrographic  Office  "Hydrographic 
chart"  April-June,  1962) 
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APPENDIX 


Since  many  examples  for  dynamics  of  the  meanders  in  the  text 
are  taken  from  the  Kuroshio,  its  behavior  is  described  briefly.  Hydrographic 
surveys  of  the  Kuroshio  after  the  World  War  II  were  mostly  limited  north 
of  25° f  although  the  area  from  the  east  of  the  Philippine  to  the  north  was 
frequently  surveyed  before  19li0.  On  the  other  hand,  the  area  south  and 
east  of  Japan  between  130°E  and  l55°E  was  surveyed  at  least  four  times  a 
year  by  five  to  eight  research  vessels  Almost  synoptically  since  1  5>0  and 
thus  seasonal  changes  aid  meanders  of  the  Kuroshio  are  well  known.  The 
historical  development  of  studies  on  the  Kuroshio  in  Japan  was  summarized 
by  Uda  (196U) - 

In  order  to  show  somefeatures  of  the  Kuroshio  near  Japan,  two  charts 
of  surface  dynamic  topography  relative  to  the  1000  d-bar  surface  in  May-June, 
19^1  and  Tray,  1955  are  reproduced  in  fig. 2  from  Ichiye's  paper  (1957). 

In  this  figure  bottom  topography  is  also  shown.  It  is  seen  immediately 
that  the  current  pattern  showed  two  different  t*rpes  to  the  south  of  Japan 
between  136°  and  lU0° •  one  type  with  a  large  cyclonic  detour  around  an 
isolated  cold  water  as  in  1955  and  the  other  trpe  with  the  flow  close  to  the 
coast  as  in  1951*  The  separation  of  tie  Kuroshio  from  the  coast  occurs  at 
around  36°N  where  the  general  coastline  of  Japan  changes  its  direction  from 
ENE  to  almost  due  north.  Immediately  after  leaving  the  coast,  the  Kuroshio 
usually  makes  a  weak  anticylonic  detour  to  the  north  between  1^2°  and  lh6°E 
as  a  first  half  wave  length  of  meanders.  The  sea  to  the  east  of  l50°E  was 
surveyed  only  occasionally  and  multiple  currents  seem  to  occur  as  in  the  Gulf 
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The  southern  detour  of  the  current  is  peculiar  to  the  Kuroshio  and 
has  no  counterpart  in  the  Gulf  Stream.  It  lasted  for  several  years,  although 
the  path  of  the  flow  changed  in  some  degrees  d '.ring  its  occurrence.  There 
are  many  specula 'ions  on  the  cause  of  this  detour  but  none  of  them  is 
conclusive.  However,  since  the  detour  occurs  just  upstream  of  the 
Bonin  Ridge  which  also  has  no  counter  part  in  the  Gulf  Stream  area,  this 
ridge  may  be  a  main  cause  for  the  detour.  In  Fig.1‘3  the  latitudes  of 
the  southern  most  limit  of  the  Kuroshio  between  135°  and  ll(.0°  were  determined 
from  temperature  distribution  at  100m  and  dynamic  topography  as  an 
indication  of  the  extent  of  the  southern  detour  and  the  result  is  plotted 
with  comments  on  existence  of  a  cold,  cyclonic  vortex  (Ichiye,  1957)* 

Also  in  this  figure  were  plotted  the  northern  limit  of  the  Kuroshio  along 
lhh°E  determined  from  temperature  profiles  and  the  schematic  trend  of  the 
water  temperature  measured  at  t.iree  coastal  stations  along  the  northern 
coast  of  Japan  (Kawai,  1955)  •  It  appears  that  the  southern  limit  changed 
with  occurrence  of  the  cyclonic  eddy  x^jith  period  of  eight  to  ten  years  while 
the  northern  limit  fluctuated  with  much  shorter  period  of  four  to  five  years. 
Therefore  in  some  periods  of  years  the  southern  and  northern  limits  were 
much  apart  while  in  others  they  were  close  together.  A  qualitative 
explanation  of  such  coupling  and  decoupling  process  was  discussed  by 
Ichiye  (1957)*  Moriyasu  (1963)  published  series  of  charts  showing  the  axis 
of  the  Kuroshio  between  131°  and  ll|0°S  and  cold  water  mass  accompanying 
Ttfith  its  cyclonic  detour  from  1931  to  1961. 


. 

. 

. 

. 

r  : 

.  . 


APPENDIX  -  3 

Masuzawa  (1^61)  determined  the  oath  of  the  axis  of  the  Kuroshio 
between  130°and  150°E  from  surface  velocity  measured  with  GlK  and  from 
subsurface  temperature  measured  with  BT  from  195h  to  1959*  His  results 
are  shown  in  Eig.llu  Except  in  1959,  the  oath  west  of  136°e  was  almost 
constantly  close  to  the  coast  and  did  not  show  much  fluctuation.  The 
southward  detour  between  136°  and  liiO0  started  in  the  middle  of  1953, 
disappeared  in  1956  and  reappeared  in  1959«  The  meanders  east  of  l!i5°E 
fluctuated  more  rapidly  and  particularly  a  cyclonic  southward  bend  had 
much  stronger  curvature,  large  amplitude,  and  shorter  duration  than  the 
one  between  136  and  lliO0,  suggesting  that  t  :ese  meanders  may  be  generated 
by  an  instability  of  a  free  jet.  The  northern  limit  of  the  Kuroshio  between 
lU2°  and  lItO°E  and  the  southern  limit  between  136°  and  139°E  are  ■'lotted 
also  by  Masusawa  (lf6l).  He  showed  that  shift  of  the  northern  and  southern 
limits  was  almost  in  the  same  direction  until  1959,  though  the  northward 
shift  o"  the  northern  limit  preceded  the  shift  of  the  southern  limit  by 
about  a  year.  The  occurrence  of  the  southern  detour  in  1959  was  discussed 
by  Moriyasu  (1961)  as  four  successive  stage  based  on  temperature  at  100m 
and  GEK  data  collected  by  different  cruises  between  May  10  and  June  9»  In 
the  first  stage,  May  10  to  II4.,  the  cold  eddy  of  about  200  km  length  in  the 
east-west  direction  and  100  km  width  in  the  north-south  direction  with 
its  center  at  31°N  and  133«5°E  was  first  detected  b  •  its  18°C  isotherm 
at  100  m  depth,  where  the  main  current  of  the  Kuroshio  flowed  along  the 
20°  isotherm.  In  the  second  stage.  May  17  to  20,  this  eddy  moved  eastward 
and  its  center  was  located  at  32.5°N  135°E.  In  the  third  stage,  May  22  to  26 
the  detour  of  the  current  between  13h°E  and  137°E  was  almost  the  same  as  in 
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the  second  stage  but  the  upstream  side  vest  of  13h°E  was  found  closer  to 
the  shore  than  before.  In  the  fourth  stage  from  May  29  to  June  9  the 
center  of  the  eddy  was  located  at  33°N  137*5°E,  further  eastward  from  the 
location  in  the  third  stage.  The  center  of  the  eddy  moved  slightly  eastward 
in  August  to  13'3°E  as  indicated  in  Fig.lij. 

The  maximum  surface  velocities  measured  with  GEK  are  averaged  over  a 
period  between  19 and  1 959  by  Masuzawa  (1961)  and  between  195>5>  and  1963 
by  Uda  (1961)  for  segments  with  different  longitudes  between  130°E  and 
15>5°E.  These  averged  values  are  from  2.6  to  3*7  knots  according  to  the 
former,  while  TJda's  values  are  generally  lower  by  0.2  to  O.h  knots.  But 
both  authors’  results  indicated  that  the  averaged  maximum  speed  is  larger 
between  13^°E  and  133°E  and  between  llt3°E  and  lli6°E  and  smaller  between 
130°E  and  135°E  and  between  lh9°E  and  1?2°E*  This  might  be  related  to  the 
meanders  of  semi-permanent  type,  although  their  existence  is  masked  by 
those  with  shorter  neriods.  As  for  the  vertical  structure  of  a  velocity 
field  in  the  Kuroshio,  many  authors  disagree  in  their  selection  of  the 
level  of  no  motion.  The  consensus  is  that  the  current  at  1000m  still 
flows  in  the  same  direction  as  the  surface  current  with  speeds  9w  to 
30c?  of  the  surface  velocity  for  the  deaths  of  no  motion  of  1200m  and 
3000m,  respectively.  (The  1200m  depth  was  determined  by  Ichiye  (i960)  from 
comparison  of  GEK  measurements  with  the  dynamic  calculation  and  the  3000m 
death  was  obtained  by  Nan’niti  (1962)  with  direct  current  measurements 
using  neutrally  buoyant  floats.) 
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Another  interesting  feature  of  the  Kuroshio  which  has  no  counterpart 
in  the  Gulf  Stream  is  its  branch  current  in  the  Japan  Sea,  usually  called 
Tsushima  Current.  The  temperature  distribution  at  100m  depth  and  surface 
currents  measured  with  GEK  during  May- June,  1962  are  shown  in  Fig.l5* 

Hi  is  current  is  supposed  to  be  separated  from  the  Kuroshio  at  the  latitude 
of  30°N  and  between  127°  and  130°E  but  it  is  not  distinguished  as  a  strong 
boundary  current  up  to  the  nort  east  of  the  Tsushima  Straits,  owing  to  the 
shallow  depth  and  complicated  coastlines.  Inside  the  Japan  Sea,  the  main 
part  of  the  current  flows  close  to  the  coast  with  a  maximum  speed  of  1  to 
1.5  knots  and  with  a  width  of  50  to  100km.  It  grox^rs  rather  stronger  as  it 
flows  northeastwards  and  at  the  latitude  of  lil°N  the  greater  part  of  the 
current  flows  out  to  the  Pacific  side  through  the  Tsugaru  Straits.  The 
remaining  part  of  the  current  flows  further  north  and  finally  it  is  dissipated 
by  mixing  with  the  cold  water.  Between  35°N  and  l£)0N  the  current  usually 
produces  two  anticyclonic  vortices  of  dimensions  of  the  order  of  a  few 
hundred  kilometers  to  the  left  hand  side  (offshore).  On  some  occasions 
the  main  flow'  detours  clockwise  around  one  or  two  of  these  vortices  and 
thus  only  a  fraction  of  the  total  transport  moves  close  to  the  coast.  It 
is  unknown  whether  such  a  detour  is  seasonal  or  secular  or  how  these  two 
eddies  change  their  extents  and  locations.  Hydrographic  data  indicate  that 
the  current  is  geostrophic  but  its  depth  is  rather  limited  to  about  300 
meters.  Since  this  current  has  its  coastal  boundary  to  the  right  hand  side 
the  stabilizing  effect  of  the  coast  is  not  the  same  as  in  the  Gulf  Stream 
and  the  Kuroshio.  Also  its  thermohaline  structure  is  somewhat  different  from 
the  latter  currents,  since  a  permanent  thermo-cline  of  the  Tsushima  Current 
is  located  between  100  and  ISO  meters  and  the  temperature  below  300  meters  is 
almost  uniformly  0  to  1°  centigrade. 
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